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Summary 


The observational data on the interstellar polarization of starlight is shown 
to be generally favourable to a helically twisted model of the spiral arm 
magnetic field. ‘The data is, in some respects, seriously opposed to the 
normal cylindrical model. There is also evidence that a halo field dips through 
the galactic plane in interarm regions. 





1. Introduction.—It was noted by Hiltner, several years ago, that interstellar 
polarization effects tend to attain a maximum intensity for stars in galactic longi- 
tudes close to 102°. The directions of polarization in such longitudes are close 
to being parallel to the plane of the Galaxy. These facts have been interpreted in 
accordance with the theory of Davis and Greenstein, which will first be briefly 
discussed. 

According to this theory, the interstellar grains are needle-shaped. They are 
subject to hysteresis effects that cause their long axes to tend to become perpen- 
dicular to the magnetic field. They are otherwise free to rotate; in particular 
they rotate about short axes parallel to the magnetic field. Consider, now, 
initially unpolarized light.travelling toward the Earth. The light has two inde- 
pendent directions of linear polarization, both perpendicular to the line of sight. 
At any point of the path of the light we choose one of these directions to lie in 
the plane determined by the line of sight and by the direction of the magnetic 
field at that point, the other component being taken perpendicular to this plane. 
Initially, the electric field components of the unpolarized light are closely equal. 
Scattering by the interstellar grains can produce measurable inequality, however, 
for the reason that the grains do not possess equal average projections in the two 
directions—the average projection perpendicular to the plane being the greater 
by a factor (cos @)~1, where @ is the angle between the magnetic field and the line 
of sight at the point in question. This causes the electric field component per- 
pendicular to the plane to be more weakened by scattering than the component 
in the plane. 

These considerations, relating to a particular point on the path of the light, 
have application to the case where the main scattering occurs within a particular 
localized interstellar cloud, provided the magnetic field inside the cloud can be 
taken to be approximately uniform. Then we have effectively the same plane 
(defined by the line of sight and the magnetic field) at all points of the light path 
within the cloud—and if scattering is unimportant outside the cloud, other parts 
of the light path need no discussion. Scattering causes a greater weakening of 
the electric field component perpendicular to the plane than of the component in 
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the plane, except in the special case where 0=0. Write A for the latter com- 
ponent, as measured by the terrestrial observer, and B for the former. The 
degree of polarization P is then defined by 
A-B 

P= A+B’ (1) 
Now we can write A=Jexp (—acos@), B=Iexp (—a), where J is the initial 
unscattered intensity, and a is a constant for the particular cloud in question, 
depending on the size distribution and number of the grains that it happens to 
contain. Evidently P does not depend on J. Furthermore, since the observed 
values of P are <1 we can infer that in general a <1, since 6 cannot be close to 
zero in all cases. With a <1 we evidently have 


Pz ta(1—cos 6). (2) 


By determining observationally the maximum and minimum values of the 
electric field, the quantities A, B are obtained. The direction in which the 
maximum is found determines the plane of vibration containing the line of sight 
and the magnetic field. Although both the value of P and the plane of vibration 
are thus observable, the value of @ cannot unfortunately be found directly, since 
the value of a is not known—and indeed a presumably varies from one cloud to 
another. The line in which the plane of vibration intersects the sky can readily 
be determined, however. ‘This gives the projection of the magnetic field on to 
a plane perpendicular to the observer’s line of sight. A convenient way of 
representing the observational data in any particular case is to draw a segment 
of a line with length proportional to the value of Pon a map of the sky (cf. 1), the 
centre of the segment being at the position of the observed star and the direction 
being that of the projected magnetic field. 

When represented in this way, the facts mentioned in the first paragraph show 
up as segments of maximum length in galactic longitude /~ 102°, the directions 
of the segments being nearly parallel to the plane of the Galaxy. If we make 
the reasonable postulate that the values of a for the interstellar clouds possess 
no statistical correlation with respect to heliocentric galactic longitude, then 
maximum values of P in /~ 102° imply that cos@~o for this longitude. This 
would be the case if the solar system lies close to or within a spiral arm in which 
the magnetic lines of force are everywhere parallel to direction /~ 12°. 

The situation described in the previous paragraphs constitutes the main 
evidence for the frequently quoted view that magnetic lines of force within a 
spiral arm run everywhere parallel to the arm. It will be the object of this 
paper to examine how far such a view is really a necessary inference from the 
data. It will be shown that the inference is not necessary. Indeed it leads to 
an immediate and serious difficulty. 

2. A difficulty in the usual model of the spiral arm magnetic field.—If the 
magnetic lines of force are parallel to the spiral arm and if they are directed 
toward /~ 12°, then the spiral arm must itself be directed toward /~12°. This 
differs by some forty degrees from the direction indicated by 21cm hydrogen 
line-data, which suggests that the arm must be directed toward /~ 54°. 

In an attempt to evade this difficulty it has been argued that the directivity 
toward /~ 12° represents simply a local deflection of the arm, and that the 21 cm 
data refers to a larger scale directivity. ‘The polarization data itself allows this 
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suggestion to be put to test. The test does not support the notion of a local 
deflection. 

Remembering that the solar system lies on the inside of the Orion arm, the 
idea of local deflection would require directions such as /= 50° to pass immedi- 
ately through the arm. Polarization of somewhat more than half the amount at 
1= 102° would thus be expected at /~50°, the direction of polarization being 
still just as parallel to the plane of the Galaxy as it is at/=102°. But if, on the 
other hand, the arm is directed toward /= 54°, a direction such as = 50° would 
pass into an interarm region, in which case the regular behaviour of the polari- 
zations might be expected to disappear. Reference to the data for regions of 
the sky within 5° of the galactic plane (1) shows that parallelism to the galactic 
plane disappears abruptly at longitudes less than about 54°, just as would be 
expected on the basis of the 21cm evidence. It is difficult to see how the idea 
of a local deflection can be maintained in view of this apparently decisive test. 

It seems then as if either the usual model of the magnetic field parallel to the 
spiral arm or the Davis-Greenstein theory must be abandoned. Of these, the 
first alternative is much to be preferred, especially since quite independent con- 
siderations cast doubt on the usual model. In a previous paper (2) it was 
suggested that the magnetic field of the spiral arm has a helically twisted structure. 
By adopting such a picture several difficulties were avoided, particularly the 
difficulty of preserving the arms as coherent structures. On the usual model of 
a simple cylindrical magnetic field the arms are easily disrupted by the tidal 
action of the gravitational field of the Galaxy. Disruption is effectively resisted 
by a sufficiently tightly wound helical structure. 

The immediate question arises concerning a helical magnetic model, as to 
whether such a model is at all consistent with the polarization data, and in 
particular whether the difficulty discussed above can be resolved by a helical 
model. This question will be examined in the remaining sections of the paper. 

3. An extreme helical model.—The case of lines of force wound in tight helices 
will be treated. Such a case is evidently the extreme opposite to the simple 
cylindrical picture. At first sight one would expect such a radical change of 
model to lead to grave discordance with the polarization data. That this does 
not turn out to be so demonstrates very clearly the incorrectness of the often 
quoted notion that the polarization data supports the cylindrical picture. Indeed 
the extreme helical model appears to fit the data better than the cylindrical model. 

Before proceeding, it may be useful to remark that the field of the arms 
becomes twisted if there is a steady drift of material outwards from the nucleus 
of the Galaxy and if the field is highly twisted by rapid rotation in the nucleus. 
This was the picture arrived at in the earlier paper (2). 

It is worth emphasizing that the degree of twisting of the spiral arm field 
could hardly be sufficient to lead to the extreme helical model treated here if the 
velocity curve given in a previous communication (3) was adopted. For this 
velocity curve, of the order of ~ 10-100 twists would be distributed over a spiral 
arm of length ~30 kpc and this would give rise to a loosely wound helical structure. 
An extreme helical model is suggested by the velocity curve of (4) where the 
number of twists is much larger. Preliminary results of the Dutch astronomers 
have confirmed the general shape of the velocity curve of (4); however the 
rotational velocity at the inner nuclear peak may be even higher than the value of 
200kmsec~! given in (4). The number of twists in this case would be closer 
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to 1,000 than 100. This point is of importance, since it is evident that a pro- 
nounced difference in the pattern of the polarization would result between a 
loosely wound and a tightly wound helical model for the spiral arms. 

If the arms had circular cross-section, the lines of force would be expected to 
have circular turns. Because the arms are flattened in the direction perpen- 
dicular to the galactic plane, the turns will be correspondingly flattened, however. 
The turns will here be taken as ellipses with major axes in the galactic plane, and 
with centres lying on a straight line directed toward galactic longitude /= 54°. 
Since this line forms the axis of the arm, the present model agrees at the outset 
with the 21cm data. No serious discrepancy of the sort described in the previous 
section for the cylindrical model therefore arises. 





50 T T T T T T T + T T T T or 
° 

+40- 4°40 

+30 +30 


LAT'TUDE 





GALACTIC 








-20b 4-20 
-30F 4-30 
-40- 4-40 
-5 L 1 L | 1 | n 1 ! 1 1 tI _go: 





330 320 310 300 290 200 270 260 250 240 230 220 210° 
Fic. 4.—Polarization of all stars observed between longitudes 210° and 330°. 


The crucial feature of the present model is that the elliptical turns of the field 
are taken to lie essentially in planes that make an angle of ~ 40° with the axis of 
the arm. Thus the direction /= 12° is essentially parallel to the turns. Sucha 
feature makes good sense, because the gas on the inner side of an arm must tend 
to rotate around the Galaxy in a shorter period than the gas on the outer side. 
This tendency is resisted however by the helical magnetic structure but only after 
the turns have become distorted in the present manner. That such a distortion can 
prevent any relative streaming of the gas between inside and outside was shown 
in the previous paper (2). 

The polarization produced by discrete interstellar clouds located on the 
elliptical turns have been worked out for three separate cases: 

(1) Semi-major axis 300 parsecs, semi-minor axis 75 parsecs. 

(2) Semi-major axis 250 parsecs, semi-minor axis 100 parsecs. 

(3) Semi-major axis 125 parsecs, semi-minor axis 50 parsecs. 

Cases (1) and (2) were chosen to show the effect of changing ellipticity, and (2) 
and (3) to show the effect of varying scale withinthe arm. In each case the solar 
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system was taken to lie at a distance of 400 parsecs from the axis of the arm, and 
at a height of 25 parsecs above the galactic plane (i.e. north of the plane). The 
latter value is believed to be correct to within about 5 parsecs, a circumstance of 
importance since the results are somewhat sensitive to height above the plane. 
The results are much less sensitive to distance from the axis. 
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Fic. 6.—Polarization of all stars observed between longitudes 320° and 10°. 


Calculations have been made for a wide variety of points on ellipses in different 
directions. Selecting a particular ellipse the points chosen were spaced at 
angular intervals of ~20° around it. In each case about ~ 30 ellipses were taken, 
whose centres were spaced at equal intervals of distance 50 parsecs measured 
from the plane containing the Sun and perpendicular to the arm. These were 
taken out to a distance of ~ 1,000 parsecs in both directions along the arm—i.e. 
toward /=54° and /=234°. ‘The important assumption was made that the field 
within the clouds always possessed the local directivity of the helical field—that 
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is to say, was directed essentially along the tangent to the particular ellipse under 
consideration. The effect of this assumption will be discussed in the following 
section. 

The calculations determine the direction of projection of the magnetic field 
on a plane perpendicular to the line of sight—i.e. on the sky—and also of the value 
of 1—cos@. Results are plotted in Figs. 1, 2, 3 for the respective cases. The 
effect of each cloud is represented as a segment, the direction of the segment 
being the direction of the projected magnetic field and the length of the segment 
being proportional to 1 —cos 0. 

Distances are not marked in the figures, but a distinction has been indicated 
between a cloud on the inner (near) side of an ellipse and one on the outer (far) 
side—the former are continuous lines, the latter are dotted. In this respect the 
figures resemble the observational data, no precise distances being determined 
for the observed polarization effects. 

For comparison, the essential features of the observed polarization data are 
reproduced in Figs. 4, 5, 6 (corresponding to Figs. 1, 2, 3 of the paper by Hall (1)). 

4. Discussion of results as compared with observations of polarization in the 
Orion arm.—The following points of broad agreement may first be noted: 

(1) Polarization effects are confined to /> ~60°. The effects of the Orion arm 
cease at about this longitude. The mosaics given by Hall (1) show almost no 
cases of polarization at /~ 60°, although at /~70° a large number are found. 
Longitudes < ~60° correspond to the interarm region (assuming no distances 
> ~1kpc), as already noted above. The calculations do not apply to the inter- 
arm region (see remarks of Section 6). 

(2) For ~70°</l< ~130°, polarization is strong, and statistically there are many 
more cases where the polarization is nearly parallel to the galactic plane than there 
are at a large angle. 

(3) For ~170°</l< ~220°, the polarization is statistically much weaker, and 
there is a marked tendency for the few stronger polarizations to turn at right 
angles to the galactic plane. The observations show few cases of polarization 
close to the galactic plane. At galactic latitudes |b] > ~5°, some polarizations 
are found, many indeed nearly normal to the galactic plane. 

What is not shown by the observations are the criss-cross polarizations at 
steep angles for b= —2° and for longitudes ~go°</<~140°. In part, 
obscuration may have limited the data immediately to the south of the galactic 
plane, for there is a notable absence of observed cases precisely at b> —2° for 
~105°<l< ~140°. , A perhaps more important point is that clouds probably 
do not stay at fixed positions within the arm. Very likely they oscillate up and 
down through the galactic plane. It is plausible to suppose that the field within 
a cloud stays more or less fixed within it as the oscillation proceeds, the direction 
of the field not agreeing with the local field at any particular moment but with the 
field at the place where the cloud happened to be first formed. Thus the calculations 
refer more to the initial situations at the moment of formation of clouds than to 
the ensuing state of affairs. It is easy to see that the general effectof the ensuing 
oscillation is simply to produce a mixing of polarizations with respect to galactic 
latitude, but not significantly with respect to galactic longitude. The mixing 
has the effect that clouds forming appreciably above and below the plane will 
sometimes be found near the plane. These clouds produce polarizations 
essentially parallel to the plane for ~70°</< ~140°. Indeed a general mixing 
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of polarizations from high latitudes with those at b= ~ — 2° immediately destroys 
the dominance of the criss-cross effect at steep angles referred to above. 

It can be argued, however, that essentially no criss-cross effect at all is found 
in the observed polarizations for ~go°</< ~105°, and that according to the 
present model some should be found, even though their comparative frequency 
is weakened by the cloud motions. This criticism is certainly valid on the present 
picture unless all the polarizations observed in this special region owe their 
origin to clouds formed appreciably away from the galactic plane. Perhaps over 
this particular longitude range and over the time interval of relevance, no clouds 
happen to have formed near the plane. 

Another possibility is that the polarization is not so much produced in discrete 
clouds as in a general medium distributed along the line of sight. This would 
yield an averaging effect, diminishing steep angles to the galactic plane. The 
disadvantage of any such averaging is that it makes it more difficult to understand 
the marked variations in the direction of polarization often found for stars that 
are very close together on the sky. An attractive feature of the present model 
is that it yields variations of the sort that are actually found in many directions. 

5. Polarizations toward the galactic centre.—The effect of the solar system 
being on the outside of an arm would be to turn the figures through 180° for an 
inner arm, if the position of the solar system with respect to the arm were the 
same, and if the arms were not curved. The latter effect in particular must tend 
to squeeze the whole polarization pattern into a somewhat smaller range of 
longitude for an inner arm than for an outer arm. Nevertheless it is to be expected 
that the arm immediately interior to the solar system should produce a generally 
similar polarization effect but turned through an angle of the general order of 180°. 
Strong polarization directed predominantly toward the galactic plane should 
thus be found near /~ 280°, and this in fact appears to be the case. 

A new effect arises toward the galactic centre, however. In directions near 
327° several arms may produce polarization of distant stars. The pattern of 
polarization from ~ 320° to ~ 350° is especially interesting (see Fig. 6) in that 
it does show a strong criss-cross effect at steep angles to the galactic plane. This 
effect, worringly absent in the Orion arm, is indeed most strikingly displayed 
toward the galactic centre. Without some form of twisted field it seems 
impossible to explain such an effect. This is a further feature in marked 
contradiction to the simple cylindrical model of the magnetic field. 

6. Polarizations in interarm regions.—The sudden change in the polarization 
pattern at /< ~60° has already been remarked on. Polarizations become small 
and fairly randomly distributed. With the axis of the arm directed in /=54°, 
this is precisely to be expected—the interarm region is encountered for stars 
at distances < ~1 kpc at a slightly greater angle than /= 54°, simply because the 
solar system is offset some 400 parsecs fom the axis of the arm. 

For the same reason in the opposite direction we pass into the interarm region 
at an angle somewhat less than 234°, say 1< ~230°. Hence we expect the 
polarizations to change markedly also at about /= 230°, and this in fact is exactly 
the case. Once again the polarizations become smaller and randomly distributed 
(cf. Fig. 4). 

A further highly interesting point is that both at /~60° and at 1~230° 
polarizations are found in high galactic latitudes up to |b|~30°. This would be 
entirely consistent with the idea developed in the previous paper (2) that the lines 








472 The effect of a helical magnetic field on the polarization of starlight Vol. 122 


of force of the halo field cross through the galactic plane in the interarm regions. 
It seems that the effect of the halo field looping down around the arm can be 
decisively seen at /~ 65° and again at/~235°. A field at high galactic latitudes 
can also be seen toward the galactic centre, as would be expected if the halo field 
also loops around the inner arms. 

7. General conclusions.—The assumption of a highly twisted spiral arm field 
is generally consistent with the polarization data, although over a limited range 
of longitude from about go° to 110° the almost complete absence of any cases 
where the direction of polarization is inclined at a steep angle to the galactic 
plane raises an important difficulty. This could be due to the simple and rather 
extreme nature of the model, or it could be due to the particular circumstance 
that all nearby clouds in these directions happened to have condensed at 
appreciable distances from the galactic plane. 

The polarization picture would be essentially different for a spiral arm field 
with a low number of twists—such as would result from the velocity curve of (3). 
Such a picture would clearly display the very loose helical structure of the model, 
and conflict with the optical data. 

There is strong evidence that at /< ~60° and /> ~230° the line of sight 
passes into interarm regions. 

The strong criss-crossing of the polarizations toward the galactic centre 
supports a twisted model for the spiral arm field, and contradicts the cylindrical 
model, as does the interarm situation. 

There is evidence that a field descends from high galactic latitudes through 
the interarm regions. This could well be the halo field crossing the galactic 
plane. There is similar evidence in directions toward the galactic centre. 

Acknowledgments.—I wish to express my gratitude to Professor F. Hoyle for 
his encouragement of the present work, and also to the staff of the Cambridge 
Mathematical Laboratory for providing computing facilities on EDSAC 2. 
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A NOTE ON EQUATORIAL ACCELERATION 
IN A MAGNETIC STAR 


L. Mestel 


(Received 1961 February 8) 


Summary 


A rotating star has initially a purely poloidal magnetic field. Meridional 
circulation of matter causes non-uniform rotation, which in turn generates a 
toroidal component of field and so a toroidal magnetic force. The possible 
steady states are studied; it is shown that if the circulation speed is less than 
the local Alfvén speed (defined by the poloidal field), then flow along a field 
line parallel to the surface will cause a steady increase of angular velocity 
towards the equator. 





Consider a rotating star, with a poloidal magnetic field symmetrical about 
the axis of rotation. In a steady state, without meridional circulation of matter, 
the angular velocity is very nearly uniform along a line of force (1). Now let 
there be a maintained meridional circulation of matter. The convection of 
angular momentum causes a strong variation of angular velocity along a line 
of force, and so generates a toroidal field. The resulting toroidal force changes 
the angular momentum of circulating elements. Further, the change in the 
toroidal field at any point depends not only on the non-uniform rotation along 
the poloidal field, but also on the convection of the toroidal field by the circulation. 
We now suppose that a steady state has been reached, and look for the distribution 
of angular velocity and toroidal field strength. 

We write the magnetic and velocity fields 


H=H,+H, (1) 
V=V,+Vj=V, t+ at, (2) 
where the suffices p and ¢ refer respectively to poloidal and toroidal components, 
the unit vector t is in the positive toroidal direction, Q is the angular velocity, 


and @ is the perpendicular distance from the axis. In the present note we assume 
the gas to be a simple perfect conductor, so that in a steady state 


V,(¥,H)=V,[(vp + ¥:),(H, + H,)] =0. (3) 

Thus 
V,(¥p,H,) =0; (4) 
V,.(v,,Hp + ¥p,Hi) =o, (5) 


separating out poloidal and toroidal components respectively. As the system 
is axially symmetric, the toroidal vector v,,H, cannot be the gradient of the 
(single-valued) electrostatic potential; hence from (4) 


v,=«H,, (6) 
where « is a scalar. Introducing cylindrical polar coordinates (@, ¢, 2), with 
origin at the centre of the star, and the rotation axis as z-axis, we have from (2), 








474 L. Mestel Vol. 122 


(5) and (6) 
~ (@QH,—«HH,)— = (HH, -GOH,) 


- (a- ey | 5 (@H)+ = (oH,) | 


ill 3 ah. 
+0 (4.35 +H,z.) (o- oe) << 


V.H= aE (aH,)+ 2 (oH) | “, (8) 
(7) reduces to 


so that 
QQ— — =a, (9) 
where « is constant on a particular field line. When «=o (no motion), (9) reduces 


to Ferraro’s law Q=a. 
From the equation of continuity 


V .(pv)=V. (pxH) =H. V(px)=0, (10) 
or 
pk=, (11) 
n being constant on a field line. The toroidal equation of motion can be written 
é [A] = pv. V(Qa%); (12) 
47 t 


the couple exerted by the magnetic force on unit volume is balanced by the 
rate at which angular momentum is transported out of this volume by the 
circulation. H, could be zero only if Qa? were constant on a field line, and 
then (9) would be violated. With 


(V,H),=- +5 (aH) (13) 


— (H).= £5 (@H,), (14) 
(12) becomes, with the help of (11) 
Bs | Hoa +H =, |x )—pxH .V(Q0*)=H.V [= - ex? | =0; (15) 
4m, “da * Oz ° 47 : 
or 
- sre + prQa? = — £, (16) 
where again f is constant on a field line. 
Equation (16) is easily interpreted. The flow of matter parallel to the poloidal 
field transports angular momentum 
p0, 267A = (pxQa*)H,A (17) 
per second along a poloidal flux tube of infinitesimal area A. This quantity 
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is not constant along the flux tube, but a steady distribution of angular momentum 
is maintained through the transport by the magnetic stresses of angular momentum 
oH, 

(- =) H,A, (18) 
as shown by Liist and Schliiter (2). As (H,,A) is constant along a flux-tube, 
we arrive at (16). 

Combining (9), (11) and (16), we find 


a(:- #7) = («+ 48), (19) 
oH, o(: - ee) = (B+ 42a"). (20) 


Equations equivalent to (9), (11) and (16) have been given by Chandrasekhar (3). 

To determine the variation of Q along a field line, we need a relation between « 
and 8. This is given by the condition that the closed loops of the initial poloidal 
field remain closed loops under the toroidal distortion that generates H,. The 
field lines are given by 





di add_ dz 
“= °s - 
Integrating round a loop, we have by (20) 
rs _ [dz (H,\ _ dz (B+47nad*) 
= bab= $5 (Zt) ~ bam, Gar 7 


This fixes B/« in terms of y, p and H,. If 7=0, B=o—no motion, no toroidal 
field. 

Suppose first that 47*/p = v,?/(H,?/47p) < 1 over the whole of a loop—i.e. the 
flow speed is everywhere less than the local speed of hydromagnetic waves along 
the poloidal field. Then as dz/H,>0, (22) implies that 


B+ 42nua,” =o, (23) 
where @, is the axial distance of at least two points on the loop (there may be 
more than two if the equation z=2(@) to the loop is more than doubly-valued). 
By (20), H, changes sign at @=@,; as it has arisen by toroidal distortion of H,, 
clearly H, cannot have the same sign all the way round a loop. Because of the 
factor #* in the denominator of (23), the points ®, will be near the axis if the loop 
considered comes close to the axis. 

Substituting (23) into (19), we have 


1-47 @,\? 
Q ~p \@ 
a 


- (=) i (24) 
p 


Consider now a field line which has a segment nearly parallel to the stellar surface, 
so that p is constant along this segment. Then as long as (477?/p)(@,/@)* <1, 
(24) predicts that Q/« is positive, and that it increases with @ As 477*/p<1, 
Q/« is certainly positive for @>,; and if 477n?/p<1, we may expect Q/« positive 
along the whole of the horizontal segment, as @, is likely to be a point near the 
axis. 

















L. Mestel 
If 477?/p>1 along a whole loop, we can still define @, by (23), but (24) is 


best written 
4mm® (@.\* _ 
Q p \@ 
a = 2 , (25) 
(=) 
p 


This predicts a decrease in Q/« along a horizontal segment. |Q| decreases along 
the whole length if 477?/p>(@,,,,/@,)"; otherwise there is a point of zero 
rotation on this segment, with |Q| increasing on either side. When 477?2/p>1, 
(25) predicts Q@?=constant on a streamline, as expected. 

Finally, suppose that (1—477?/p) changes sign at two points P, and P, of 
a loop (the discussion is unaffected by the restriction to two such points). Then 
if our steady state is to be meaningful, we must have simultaneously : 
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B+4rnad2=0;  — 477?/p=1; (26) 
and so also 
nB _ 
_ Gn) 


p, being the density at @,. For if (26) and (27) do not hold, both H, and Q have 
singularities, yielding infinite rotational and magnetic energy densities at @,. 
(It is true that an infinite rotational shear is inconsistent with the implicit neglect 
of viscous forces in (12); but stellar viscosity is so small (4) that a steady state 
that included viscous forces would still result in local centrifugal and poloidal 
magnetic forces far larger than the equilibrium of the star could tolerate.) 

Equations (26) and (27) imply a density field distorted from the approximately 
spherical, for at P, and P, both #, and p, would have to be the same. However, 
we can rule out this case by substituting (26) into (22): 


af dz (#-@,2) _ 
26 ah =0, (28) 
PoP 


Thus as numerator and denominator of the integrand have the same sign, 
a= 8=o0: asteady state is possible only if the angular momentum in this particular 
flux-tube is zero. We conclude that in general a steady state is possible only if 
the circulation is either ‘‘sub-Alfvén”’ or ‘‘super-Alfvén’’ everywhere along 
a field line. 

If / is the perpendicular distance between neighbouring poloidal loops, the 
angular momentum within the flux-tube defined by these lines is 





wera 8 
p pQ&*(27al ds) = 20 (H,al)a b e m.. ds (29) 
Pp v (: uy ee) ’ 
p 
where ds is the scalar element of length. The transition to our steady state from 
the “‘initial’’ state, with no circulation and no toroidal field, involves merely a 
redistribution of angular momentum along flux-tubes; hence given v (satisfying 


(11)) we can by (22) and (29) express « and f in terms of 7, H,, and the prescribed 
angular momentum within the flux-tube. 
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As a rough numerical example, suppose 477?/p<1; then |AQ/Q| along a 
horizontal segment from @, to @, is, by (24), approximately 


(477n°/p)[(@_/@1)’ — (@-/@g)), 
which is ~ 477?/p if d.>a@,~ @,. 
Taking for v the velocity of migration of the sunspot zone (~ 2 x 10? cm/sec) 
and for p a sub-photospheric value (10~’g/cm*), |AQ/Q|~5 x 10-*/H,”. An 
equatorial acceleration of 15 per cent requires H,, to be ~0-6 gauss. 


Discussion 


The prediction of equatorial acceleration is sufficiently remarkable to justify 
this preliminary note. However, it is clear that many problems remain. We 
have taken as given a velocity field v, that satisfies (6) and (11); in fact, the 
circulation in non-convective regions is determined by the non-spherical distortion 
of the thermal field by the centrifugal and magnetic forces (§). Over the bulk 
of the star the circulation is slow, and will generate only very slight departures 
from isorotation, and very weak toroidal fields; but near the surface where the 
density is low the circulation is in general much more rapid (6), and should have 
important hydromagnetic effects. 

If we regard the field H,, as at our disposal, and exerting negligible poloidal 
forces, we may compute v, in terms of 2 and H,, and then choose H, so as to 
satisfy (6) and (11), so that steady state conditions are fulfilled. Even so, as by 
(19) and (20) Q and H, depend partly on H,, there may still result a difficult 
self-consistent problem. In a more realistic problem, the initial field H, would 
be arbitrary. It should be noted that the thermally driven circulation is 
‘‘inexorable’’, being fixed by the Q and H, fields: even with the energy 
density of the circulation less than that of H,, the circulation must cause 
progressive distortion of H,, as long as the reaction back of H, on the thermal 
field can be ignored. For the steady state theory to be relevant, the coinciding 
streamlines and poloidal field lines must be closed loops within the body of the 
star. If, for example, the initial field is of roughly dipole type, with field lines 
entering and leaving the surface, a steady state of circulation will not be reached 
until the field has been greatly distorted, and the internal and external parts 
of the field separated by Ohmic diffusion. It is not clear that this will occur 
before the distorted field H,, exerts forces strong enough to react back on the 
thermal field, and so on the circulation, which might even be reduced to zero. 

The treatment has been based on a simplified form of ‘‘Ohm’s law’’ for 
a plasma, which neglects the Ohmic, Hall and thermoelectric fields (4). As is 
well known, by Cowling’s theorem an axially symmetric system of finite 
conductivity cannot be strictly steady; near the O-type neutral points that the 
field H,, must possess, the Ohmic term always dominates over the induction 
term, so that the field must decay. However, for large-scale fields, the decay is 
slow and we may expect the Q and H, fields steadily to adjust themselves to H, 
as it changes. [In a convective region, a Bullard-type dynamo may operate: 
the inevitable fluctuations about strict axial symmetry may allow the poloidal 
field to be regenerated from the toroidal (7). Our axially symmetric treatment 
could then be regarded as a smoothed-out version of the actual non-symmetric 
processes. | 

On taking the curl of Ohm’s law, it is found that the neglected terms have 
very little effect on (9g). However, they do affect the flow of current parallel 
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to H,, so that (22) is replaced by a new (f, «, 7) relation, in a time-scale of the 
same order as the decay time of H,. Thus the initial closed loops of field will 
not remain so. In fact, the condition (22) is a rather special one; the hydro- 
magnetic and dynamical conditions (9) and (16) can be satisfied, given the velocity 
field and the angular momentum per flux-tube, by any (8, «, 7) relation (provided 
the star is not given more disruptive energy than it can bear). Thus we could 
imagine our initial field (before circulation starts) as having both poloidal and 
toroidal components, but exerting no torque; when the circulation starts, an 
extra toroidal field is generated, exerting torque. But again the field so built 
will not satisfy the poloidal component of Ohm’s law, and so it will gradually 
change its structure. Further treatment will be given in a later paper, with 
special reference to the equatorial acceleration problem, and the amount of energy 
stored in the rotational and the toroidal magnetic fields. For the moment we 
just wish to emphasize that equatorial acceleration is not tied to the existence 
of field lines closed in space. It is condition (23) that yields the acceleration 
for sub-Alfvén flow, and this by (20) merely requires that H, should change sign 
at points distant #, from the axis; this will occur in infinitely many fields 
besides that determined by (22). 

Finally, we should note that although in convective zones large-scale meridian 
circulation is superimposed on the small-scale turbulent motions (8, 9) the 
immediate application of the results on equatorial acceleration to the solar case 
is somewhat premature: the existence of the solar cycle clearly indicates that a 
periodic rather than a steady solution of the equations should be sought. We may 


conjecture that the strong sub-photospheric toroidal magnetic girdles, which I, 
are presumably responsible for sunspots, generate circulation through the 
disturbance they themselves exert on the thermal field, and are dragged by the n 
circulation. A periodic rather than a steady solution would therefore seem more ol 
plausible, but it is not clear that such a solution would yield a constant equatorial g 
acceleration. tc 
Note.—In a private communication, Sir Edward Bullard has commented u 
on the similarity between the work of this paper and that of himself and his re 
collaborators on the westward drift of the Earth’s magnetic field (10). Ce 
() 
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THE RADIO EMISSION FROM NORMAL GALAXIES 
II. A Strupy oF 20 SpiraLs aT 158 Mc/s 
R. Hanbury Brown and C. Hazard 


(Received 1961 March 3) 


Summary 


A survey has been made of the radio emission at 158 Mc/s from twenty 
bright spiral galaxies. A radio index R has been defined by the equation 
R=m,— pg, where m, and mpg are the radio and photographic magnitude 
respectively. The results indicate that R is a function of the inclination and 
galactic latitude of a galaxy and it is suggested that this effect is due to the 
absorption of light. ‘The mean value of R for spirals of type Sb and Sc is 
+1°3 with an r.m.s. dispersion +0°7. No significant difference was found 
between the mean values of R for the two types of galaxy. 





1. Introduction 

A survey of the radio emission from some bright normal* galaxies has been 
made at 158 Mc/s using the 250 ft steerable paraboloid at Jodrell Bank. The 
object of this survey was to establish the integrated radio magnitudes of these 
galaxies, and in particular to investigate how the ratio of the total light to the 
total radio flux varies among galaxies of the same type and also how it depends 
upon type. In Part I (1) we described the apparatus used in this survey and 
reported the results obtained on M31 and M33; the present paper (Part II) 
contains the results on twenty spiral galaxies of types Sb and Sc. The 
observations of other types of galaxy will be published in a further paper 
(Part II1) together with a discussion of all the data. 


2. Method of observation 

A description of the apparatus was given in Part I together with the method 
of observation of M31 and M33. For the remainder of the galaxies the 
observations were carried out as follows. First, a preliminary survey was made 
in an effort to detect radio emission from a particular galaxy. In this survey 
two sets of scans were made through the expected position, one set of scans being 
made at constant right ascension and the other at constant declination; each 
scan was repeated several times as a check on its reliability. In some regions 
these scans sufficed to give the position and intensity of the radio source 
associated with the galaxy, and an example of such a region is illustrated in Fig. 1 
which shows two scans made over NGC 253. In other more confused regions 
it was found necessary to survey a considerable area around the expected position 
as an aid to the interpretation of the observations. 


* In this context a normal galaxy, following Baade and Minkowski (2), is taken to be a galaxy 
which appears to,be normal when examined photographically. 
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Fic. 1.—Two scans in declination across NGC 253. These scans were both made at a constant 
right ascension (ooh 45™-1) with a beamwidth of 2° at 158 Mc/s (uly, 1958). 


Before each set of observations the equipment was calibrated on the source 
I.A.U. 14N5A and all the magnitudes in the present paper were measured 
relative to that source. The scale of radio magnitude (m,) is defined by the 
equation 

m, = — 53°45 — 2°5 log Sisg (z) 
where S,;,w. m~*(c/s)~ is the flux density of a source at 158 Mc/s. Following 
Part I we have taken the magnitude of 14N5A to be + 6-86. 


3. Results 

Observations were made on the twenty galaxies listed in Table I. All these 
galaxies are brighter than m,,= 10-6 on the scale of Holmberg (3) and were in 
convenient positions to be surveyed during the autumn of 1958 when this work 
was carried out. 

The results are summarized in Table I, which gives the positions and apparent 
magnitudes of the radio sources; some detailed remarks on the individual 
galaxies are contained in the notes. In fifteen cases a radio source was found 
close to the position of the galaxy, and the measurements are sufficiently precise 
to allow an identification to be made. It is of course possible that some of the 
suggested identifications represent chance coincidences with discrete sources; 
but at 158 Mc/s there is about one source per 20 square degrees with an intensity 
above the limiting sensitivity of the present survey and it is therefore unlikely 
that more than two of the identifications suggested in Table I represent a chance 
coincidence. In the remaining five cases the isophotes of the region around the 
galaxies were too complicated to be interpreted unambiguously and a definite 
conclusion would require superior resolving power; in four of these cases there 
was evidence of a radio source close to the position of the galaxy but the position 
could not be established with satisfactory precision. When no radio source 
could be identified with a galaxy an upper limit to the radio emission was 
estimated from the records and this limit is given in Table I. 


Tasrie I 


Radio data at 158 Me|s 
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Notes on Table I 
NGC 224. See Part I for a detailed discussion. 
NGC 253. The radio source is clearly defined. It has also been reported by Mills (9). 
NGC 598. See Part I. 


IC 342. Although the radio source lies on a steep gradient in both right ascension and 
declination it is clearly defined. 
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NGC 2403. This galaxy lies in a difficult region to interpret. On right ascension scans 
it lies in a sharp minimum of the background radiation and the radio source cannot be 
observed distinctly. The records indicate the presence of a source in the correct position. 
The declination scans are confused by nearby sources. 


NGC 2683. Region complicated by nearby sources. 


NGC 3031. ‘The source lies in a region of enhanced emission extending about fifteen 
degrees in both declination and right ascension which contains a number of weak sources. 
Observations are confused by a weak source about half a degree to the north of the galaxy. 
This confusing source has been shown by Palmer (6) to have a diameter of about 2 minutes 
of arc. The source in the position of NGC 3031 appears to have a diameter greater than 
40 minutes of arc and it seems reasonable to assume that this source may be identified with 
NGC 3031. Its intensity has been estimated by subtracting the calculated contribution of 
the small diameter source from the flux density observed, an allowance being made for the 
different positions of the sources in the beam; the intensities in each case were measured 
relative to the source I.A.U. 14N5A in order to avoid errors in calibration. 


NGC 4236. This region is complicated by nearby sources. 


NGC 4244. Declination scans show a clearly defined radio source close to the position 
of the galaxy which is situated on the edge of a localized region of enhanced intensity. 
The right ascension scans were of poor quality and it was not possible to measure an accurate 
right ascension. ‘The identification is therefore rather uncertain. 


NGC 4258. ‘The radio source was well resolved on right ascension scans but the 
declination scans are too confused to permit the determination of an accurate declination. 


NGC 4490. ‘The radio source is well resolved in declination but the measurement of its 
right ascension is complicated by the region of enhanced intensity in which it lies. 


NGC 4565. A well defined source on a moderately steep gradient in right ascension. 
NGC 4631. <A well defined radio source. 


NGC 4736. The radio source is well resolved in declination but in right ascension the 
records are confused and it was therefore not possible to make accurate measurements of 
this coordinate. A provisional identification. 


NGC 4826. Although near a minimum in the background radiation this is a clearly 
defined radio source. 


NGC 5055. This region is complicated by nearby sources and gradients in background 
radiation. ‘There appears to be a radio source at the correct declination but it is to the 
west of the visible nebula. However, the interpretation of the records is subject to con- 
siderable uncertainty and there is some evidence that the observed source is a double source, 
the two components having approximately the same declination but differing in right 
ascension. One of these components may be due to the galaxy. 


NGC 5194. Although the observations are complicated by the position being near the 
centre of a region of enhanced intensity the radio source is well defined. 

NGC 5236. ‘This region is confused by a nearby source and a steep gradient, but the 
radio source in the position of the visible nebula can be clearly observed. 

NGC 5457. ‘This region is complicated by the presence of the intense source 14Ns5A. 
In the present survey the radio source associated with the galaxy and 14Ns5A are well 
resolved. 

NGC 6946. This galaxy lies in a pronounced minimum of the background radiation 
which accounts for the uncertainty in the measured intensity. 
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4. The integrated radio magnitude 

The radio magnitude m,(@), tabulated in Table I, was calculated from the 
observed peak deflection. It is therefore a function of the aerial beamwidth @ 
and the angular size of the source. The integrated radio magnitude m,, which 
is also shown in Table I, is a measure of the total flux received from a source and 
is independent of beamwidth. It is obtained by integrating the flux over the disk 
of the source, and therefore corresponds to the magnitude which would be 
observed with a beamwidth large compared with the angular size of the source. 

In the case of M31 and M33 the values of m, were found by integrating the 
detailed isophotes of the region around each source as discussed in Part I. For 
the remaining sources in the table, which have a smaller angular size than the beam, 
the differences between the observed magnitudes m,(@) and the integrated 
magnitudes m, are small, and we have calculated these differences by estimating 
the angular size of the radio sources as follows. 

It was shown in Part I that in both M31 and M33 the majority of the radio 
emission arises in an extensive corona which envelopes the visible nebula and 
we have assumed that other galaxies of types Sb and Sc possess similar coronae. 
This assumption is supported by an examination of some records made with an 
interferometer at 158 Mc/s (6). ‘These records cover the regions around 
NGC 3031, 4736 and 5194 and were made with baselines of 60A and 330A. 
They show that the majority of the radio emission from both NGC 4736 and 
5194 arises in a source with a diameter of about 25’, while in NGC 3031 the 
corresponding diameter exceeds 40’. These diameters have been estimated on 
the assumption that the sources are circular disks with a gaussian radial variation 
of brightness. Although these results are approximate they indicate that for these 
three galaxies the ratios of optical to radio diameter are similar to those found 
for M31 and M33 and they support the assumption that all Sb and Sc galaxies 
have extensive coronae. 

In order to estimate the angular sizes of the sources shown in Table I we 
have therefore assumed that every galaxy has a corona whose ellipticity and size 
relative to the visible nebula are similar to those found in M31 and M33. 
Following the results in Part I we have taken this corona to be ellipsoidal with 
an axial ratio of o-5 and have assumed that the principal planes of the corona 
and the visible nebula coincide; the apparent angular size of the major axis 
of the corona between half intensity points has been taken as 1-3 times the 
maximum extent of the visible nebula as measured by Holmberg (3). On these 
assumptions the angular size of the radio source associated with any particular 
galaxy can be estimated from the angular size of the visible nebula. The major 
axis is simply proportional to the maximum extent of the nebula as stated in the 
previous paragraph, while the minor axis can be calculated from the inclination. 
The inclination may be calculated to a sufficient degree of accuracy on the 
assumption that the visible nebula is also spheroidal, from the equation 

sin?i = ac ae me ; (2) 

I—p 
where i is the angle of inclination of the principal plane of the visible nebula to 
the line of sight, p is the axial ratio of the spheroid and 6/a is the axial ratio of 
the nebula image. We have calculated the inclination of all the spirals in Table I 
assuming p=0'2 (3); in cases where (b/a)<o-2 we have taken the inclination 
to be zero. These values of inclination have been used to estimate the apparent 
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angular dimensions a’ x b’ of the radio sources, again by means of equation (2), 
taking p=0°5 which is the value appropriate to the coronae. 

The estimated angular dimensions of each radio source have been used to 
calculate the integration correction Am/(6,b’,a’) which was applied to the 
observed magnitude m,(@) to find the integrated magnitude m,. It was assumed 
in this calculation that the distribution of intensity across the elliptical disk of 
a radio source follows a gaussian variation along both axes, and that the aerial 
beam is circularly symmetrical with a gaussian profile. On these assumptions 
it is simple to show that 

Am(8, b’, a’) =m, —m,(0) = — 2-5 log [(1 +s?)(1 + s?r?) ]#? (3) 
where s is the ratio of the width of the major axis of the radio source to the width 
of the beam, and r is the apparent axial ratio (b’/a’) of the radio source. 

The calculated values of Am(6, a’, 6’) are given in Table I together with the 
corresponding integrated magnitudes m,; for M31 and M33 the integrated 
magnitudes are those measured directly from the isophotes (see Part 1). The 
table shows that, with the exception of these two galaxies, the integration 
corrections are small and in most cases less than the probable errors in 
measurement. It is therefore unlikely that the values of m, are significantly 
in error due to the finite angular size of the radio sources, serious errors could 
only be introduced if the ratio of the size of the corona to the visible nebula is 
much greater in some of the galaxies than in M31 and M33. 





5. Comparison of the radio and photographic magnitudes 

5.1. Introductory remarks.—Since we are concerned in the present work 
to establish the ratio of the total light to the total radio flux emitted by a galaxy 
it is convenient to define a measure of this ratio. By analogy with the concept 
of colour index we propose to define the radio index R of a galaxy as the difference 
between the integrated radio and photographic magnitudes. Thus 


R=m,—mMy, (4) 
where m, and m,, are the radio and photographic magnitudes respectively. 

It is clear that the value of R, as observed for any particular galaxy, will 
depend upon the absorption of light and also of radio waves. In an effort to arrive 
at values of R which can be usefully compared, we have therefore corrected the 
observed results for the effects of absorption. We have assumed, following an 
earlier analysis (7), that for all directions of observation the radio absorption 
both in our Galaxy and in external galaxies is negligibly small, and we have also 
assumed that the radio emission from a galaxy is isotropic in direction. On the 
other hand the absorption of light cannot be neglected and may take place both 
in our Galaxy and in the galaxy under observation; these two cases will be 
considered in turn. 

5.2. The effect of the internal absorption of light in a galaxy on the radio index.— 
The internal absorption of light in a spiral is mainly due to obscuring clouds 
concentrated in the plane of the galaxy and will therefore depend on its orientation 
with respect to the observer. This effect will cause the observed radio index to 
depend on the orientation of the galaxy, and hence there should be a correlation 
between the measured values of R and the apparent axial ratio b/a of the 
photographic image; if the radio absorption is negligible, as assumed above, 
then spirals seen edge-on should exhibit low values of R. We have analysed our 
results to test for this effect as follows. 








rial 
ons 


(3) 
ith 


the 
ted 
‘he 
ion 


tly 
ild 
is 


xy 
pt 
ice 


4) 


rill 
ve 
he 
an 
on 
so 
he 
th 


ds 


on. 


on 
he 


ur 

















No. 6, 1961 485 


First, we have restricted our analysis to galaxies in Holmberg’s catalogue 
so that measurements of m,, shall be homogeneous. Secondly, in order to avoid 
confusion with the effects of absorption in the Galaxy, we have farther restricted 
the analysis to those galaxies which are more than 40° from the galactic plane; 
this excludes M31 and M33 and has the advantage that for all the galaxies included 
in the analysis the integration corrections are small and are therefore less likely 
to introduce error into the following discussion, 

The ten galaxies which fulfil the necessary criteria are shown in Table II, 
and in Fig. 2 the radio index for these galaxies is plotted against the apparent 
axial ratio b/a of the associated visible nebulae. Although the scatter of the 
results in Fig. 2 is large, they indicate that the radio index decreases with b/a 
and hence depends upon the inclination of the visible nebula. 


The radio emission from normal galaxies: II 


TasBie II 
The effect of the internal absorption of light on the radio index of a spiral galaxy 
Type 
Galaxy (after bla R(,B) 
Holmberg) 

NGC 

4244 Sc+ 0°16 —1'2 
4565 Sb+ 0°18 —0'5 
4631 Sc+ 0°23 -—0'3 
4258 Sb+ 0°40 +0°4 
3031 Sb- o'4I +1°6 
4490 Sc+ 0°53 —o°6 
5194 Sc— 0°67 +0°2 
4826 Sb- 0°67 +or1 
4736 Sb—- 0*90 +06 
5457 Sc— 1°00 +1°2 


b/a is the apparent axial ratio of the visible nebula; R(, 8) is the radio index (m,—mp,), 
where i and f represent the inclination and galactic latitude of the galaxy respectively. 
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Fic. 2.—The relation between the inclination of a spiral and the radio index. 
Ordinate: The radio index R (i, B). 
Abscissa: The apparent axial ratio b/a of the photographic image of the nebula. 
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These results may be compared with the measurements of the internal 
absorption in spirals made by Holmberg (3). He gives a plot of the total 
photographic absorption as a function of the apparent axial ratio for two groups 
of spirals, the first group comprises spirals of types Sc+ and Sc— and the 
second group comprises types Sb+, Sb— and Sa. His results for the two groups 
have been averaged and the variation of absorption with b/a is shown by the 
broken line in Fig. 2; the zero of this curve has been arbitrarily chosen to 
facilitate a comparison between the results. The radio observations are in 
reasonable agreement with the photographic measurements and this comparison 
therefore confirms the suggestion that the dependence of the radio index on b/a 
is due to the internal absorption of light. It is interesting to note that similar 
radio observations, if conducted on a sufficiently large number of galaxies, 
might provide a valuable technique for measuring the internal absorption of 
light. 

In view of the agreement between the present observations and those of 
Holmberg, we have used his values to apply an inclination correction Am(i) to 
the observed values of the radio index in order to correct them all to a standard 
inclination of go°. These values of Am(i) are shown in Table III and have 
been added to the observed radio index R(i, 8) to give the radio index R(go, B) 
corrected to an inclination of go°. 

However, it should be pointed out that the observations of Holmberg show 
a considerable variation in the internal absorption between individual galaxies ; 
the broken line in Fig. 2 being an average of measurements of a large number of 
objects. It is only for values of b/a less than about 0-4 that the mean value of the 
absorption becomes significantly greater than the dispersion between individual 
galaxies. The radio measurements show a similar effect, a significant change of 
radio index with b/a being apparent only for the three almost edge-on spirals. 
It follows that when dealing with individual galaxies it is only worthwhile to apply 
an inclination correction when the value is of b/a less than about 0-4. However, 
for less inclined galaxies the corrected values are of importance when considering 
the average radio index of a number of galaxies. 

5.3. The effect of the absorption of light in our Galaxy on the radio index.—The 
effect of absorption of light in our Galaxy would be to decrease the observed 
radio indices at low galactic latitudes. Only two of the galaxies observed lie 
close to the galactic plane where this absorption is likely to be serious, namely 
IC 342 and NGC 6946. The average radio index of these galaxies, corrected to 
inclination 90° is —o-2 compared with the average for the remainder of +09, 
suggesting that the absorption at latitude 11° is about 1-1 magnitude. 

The results therefore suggest that the radio index increases with galactic 
latitude although the data are obviously too scanty to establish the law of this 
variation. The correction for galactic absorption found by Hubble (8) gives 
the average increase of photographic absorption relative to the galactic pole 
as 0°25 (cosec 8 —1). Thus the absorption at latitude 11° should be 1-1 magnitude 
which is in agreement with the value deduced above from the radio data. 

5.4. The average radio index of a spiral galaxy.—The radio index R(g0, 90) 
for galaxies with m’,,<9g'1 is plotted in Fig. 3 against m’,,,, the photographic 
magnitude corrected for the absorption of light, where 


yg = My — Am(i)— Am(B). (5) 
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The distribution of results in Fig. 3 suggests that the two bright galaxies 
M31 and M33 have unusually high values of R. Since both these galaxies have 
low photographic magnitudes and comparatively large angular diameters it seems 
possible that this result may have been produced by some systematic error which 
is correlated with photographic magnitude. However a re-examination of the 
original measurements and the associated corrections with this possibility in 
mind failed to explain these high values of R, and we have therefore reached the 
conclusion that both these galaxies emit less radio flux in relation to their light 
than an average spiral. 

Fig. 3 also suggests that galaxies with m,,>+8 have low values of R. 
However for these fainter objects the distribution is seriously affected by 
observational selection introduced by the limited sensitivity of the radio 
equipment. To illustrate this latter effect the limits of detection have been 
shown in Fig. 3 as a broken line; any galaxy with an uncorrected radio index 
lying above this line could not be detected. 
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Fic. 3.—The relation between radio index and photographic magnitude for spirals. 
Ordinate: The radio index R (90, 90) corrected to inclination 90° and to the galactic pole. 
Abscissa: The photographic magnitude m' yg corrected for the internal absorption of light in the emitting 
galaxy and also in our Galaxy. 


It follows from Fig. 3 that the average radio index for those galaxies which 
were detected is likely to be significantly lower than the true value for the whole 
group. In an attempt to overcome this difficulty we have made another survey 
of two of the galaxies (NGC 2403, 5055) with equipment of higher sensitivity 
and resolving power. In this survey, which was made at 240 Mc/s, the two 
galaxies were detected successfully and the resulting values of radio index are 
shown by the crosses in Fig. 3; it has been assumed in reducing these results 
to 158 Mc/s that the spectral index of both galaxies is equal to that of the reference 
source 14N5A. If the two measurements at 240 Mc/s are included, then the 
results are complete for all the sixteen galaxies with a corrected photographic 
magnitude <g:1, and the following analysis of the average radio index refers 
only to this group. 
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The average radio index for all sixteen galaxies is + 1-3 with anr.m.s. deviation 
of +08. The probable error of measurement in the index for a single galaxy 
is about +0°2; this has been estimated from the errors in the radio magnitudes 
alone and the relatively small errors in the photographic magnitudes have been 
neglected. It is therefore likely that the observed r.m.s. deviation in the radio 
index is a true dispersion in the properties of spirals and is not due to errors of 
measurement. 

As a check on the effect of the corrections which we have made for the 
absorption of light we have also calculated the average radio index for those 
galaxies in the group which have low inclinations (b/a>0-4) and which are far 
from the galactic plane (8 > 30°). The mean value of the index for these galaxies 
is the same as that for the whole group. Thus we may conclude that at 158 Mc/s 
the radio index of type Sb and Sc spirals is + 1-3 with an r.m.s. deviation of + 0-8. 

If we now separate the spirals of type Sb and Sc the corresponding mean 
values of radio index are +1-4 and +1-2, and the associated r.m.s. deviations 
are +0°8 and +0°7 respectively. Since the probable error of measurement in 
each of these mean values is about + 0:1, it follows that the present survey shows 
no significant difference between the radio index of the two types; any difference 
which may exist must be considerably less than the dispersion in the index 
between individual galaxies of the same type. 

It should be noted that, strictly speaking, this comparison refers only to 
the ratio of the radio to light flux received from the two types of galaxy, and 
does not take into account any difference in the internal absorption of light when 
received at an inclination of go°. However, according to Holmberg (3), this 
difference is small (<o-2 magnitude) and therefore the comparison may also be 
taken as valid for the ratio of radio to light flux generated in the two types of 
galaxy. 

In an earlier analysis of the radiation from normal galaxies at 85 Mc/s, 
Mills (9) found values for Sb and Sc spirals of +0-3(+0°4) and +1°6(+0-2) 
respectively but concluded that ‘‘the difference between Sb and Sc galaxies is 
hardly significant’’. The results presented here are in agreement with his 
general conclusion, but do not show a difference of 1-3 magnitudes between the 
two types. 


6. Conclusions 


The results of the present survey indicate that at 158 Mc/s the observed 
radio index of a spiral galaxy of type Sb and Sc, which is a measure of the ratio 
of the received radio flux to light fiux, is a function of the inclination of the visible 
nebula to the line of sight and of its galactic latitude; they also suggest that these 
effects are due to the internal absorption of light in spirals and in the Galaxy 
respectively. In any comparison of the total radio and light flux emitted by a 
galaxy, it is therefore necessary to take into account these absorption effects. 

When corrected for the absorption of light the mean value of the radio index 
for spirals of type Sb and Sc is +1-3 with a dispersion of +0°8 magnitude. 
No significant difference was found between the two types and any difference 
which exists is less than the dispersion between individual galaxies of either type. 
It must be emphasized that these conclusions refer only to a frequency of 158 Mc/s 
and it is possible that the difference between the radio indices of the two types of 
galaxy may be a function of frequency. 
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FLARES ASSOCIATED WITH THE 1960 NOVEMBER EVENT 
AND THE FLARE NIMBUS PHENOMENON 


M. A. Ellison, Susan M. P. McKenna and fohn H. Reid 
(Received 1961 March 29) 


Summary 


Outstanding flares occurred on 1960 November 10, 12, 15 and 20 in 
association with an active region and complex spot group which had its central 
meridian passage on November 12. 

The last three of these flares projected into space showers of cosmic-rays 
which were registered on the earth by neutron monitors at ground level 
within 30 minutes after the flare flash, bringing up to ten the number of these 
events so recorded since 1942 February 28. Geomagnetic storms of unusual 
intensity also occurred on November 12 and 13, and again on November 15. 

The flares of November 10 and 12 were photographed under excellent 
conditions on films taken at 1-minute intervals with the Lyot Ha heliograph 
at the Royal Observatory, Cape of Good Hope. These results are analysed 
and flare light-curves are plotted. 

A new phenomenon—the flare nimbus—has been found in association 

with some Class 3 and 3 + flares recorded with the heliograph. This is a dark 
absorbing halo which begins to surround the flare some few minutes after 
the filaments have reached their maximum light intensity: its duration is 
>1-2 hours and its diameter is about 300 000 km. The characteristic pro- 
perties of the nimbus are described for the five observed cases—1959 July 16, 
1960 April 1, June 1, November 10 and 12—and possible causes are discussed. 
There is a strong presumption that the phenomenon is the optical counterpart 
of the cloud of relativistic electrons whose occurrence in the flare region has 
been postulated by Boischot and Denisse in order to account for the radio 
emission continuum of Type IV. 





1. Introduction.—The spot group in latitude 27° N which crossed the Sun’s 
central meridian on November 12.2 gave rise to a remarkable series of flares, the 
terrestrial effects of which are now being referred to as ‘‘the November 1960 
event’’. The group was of Ziirich type E and contained 18 spots when it was 
first seen on the east limb on November 6. By November 12 it had developed 
into a type F group, containing 60 spots, and on November 14 its area was 1775 
millionths of the hemisphere. 

The period of greatest flare activity occurred from two days before to three 
days after central meridian passage, a situation which was most favourable for 
the reception at the earth of the accelerated particle streams of different energies. 

The following major flares were observed: 


Date Class Co-ordinates Times U.T. Source of data 
B M E 
Nov. 10 3 29°N 29°E = 102og™ =yol27M >12h30M Cape film 


12 3+ 28°No1°W s138r5™ 13230M >14h25™ Cape film 

15 3+ 26°N 35°W  o22o7™ o2t21™ o4h27™ Tokyo, Nizamiah 
Three others may be noted: an important flare, detected by its radio emission 
but missed by the optical flare patrol, which occurred on November 11 at about 
034 15™ U.T., one of Class 2 or 2+ which was seen at Tokyo on November 14, 
beginning 02" 46™, ending 05" 20™ and a third which occurred when the region 
was behind the west limb at about 20" 20™ on November 20 (see later). 
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The flares of November 12 and 15 produced showers of primary cosmic-ray 
particles which were of sufficiently high energy for their secondary effects to be 
recorded at ground level. Thus, for example, the neutron monitor at Leeds 
gave an excess counting rate of ~ 100 per cent above pre-flare level following the 
3+ flare on the 12th, and an increase of ~ 45 per cent following the 3+ flare of 
the 15th (privately communicated by Dr P. L. Marsden). 

The particles responsible for the first of these cosmic-ray events began to 
arrive at the Earth at 13" 28™, 13 minutes after the start of the flare and just before 
the time (13"48™ on November 12) of the sudden commencement (S.C.) of a 
great magnetic storm. The subsequent violent fluctuations of the geomagnetic 
field, by temporarily lowering the Stormer cut-off rigidity, appear to have allowed 
unusually great and highly irregular penetration of the cosmic-ray protons into 
lower latitudes. Another interesting feature of this magnetic storm was the 
occurrence during the recovery phase of a second sudden commencement(10" 23™ 
on November 13) of great amplitude which prolonged the storm for a further 16 
hours. It seems likely that the first S.C. may have been caused by storm particles 
which were accelerated from the flare (unobserved) on November 11 (delay time 
34-5) and the second by storm particles from the 3+ flare on November 12 
(delay time 22"-9). There followed a magnetically quiet period on the 15th 
during which the second 3+ flare began at 02"07™, and the resulting cosmic- 
ray increase started at 02" 36™ with its maximum about 1 hour afterwards. Later 
on the 15th, a third magnetic storm began with a sudden commencement at 
132 04™, 

This same active region was responsible for a third cosmic-ray event at ground 
level, which began on November 20 at 20"55™+10™ (1). The region was then 
in longitude 111° W and therefore behind the limb, but elevated portions of the 
flare and of its associated spray prominences were recorded on the Climax 
films (2). 

The magnetic storm particles, having energies of <1 Mev, arrive at the 
Earth with a time delay of >18 hours. Flare accelerated protons in the energy 
range 100-400 Mev are now known to reach the Earth following a few hours 
after flares of importance 2+ or greater, which are also accompanied by radio 
emission of Type IV (see Section 5) (3). Some forty examples of this pheno- 
menon have been studied by means of satellites (4), balloons (5) and ionospheric 
measurements of polar cap absorption (P.C.A.) (6). The protons arriving after 
the flares of November 12, 15 and 20, however, extended into a still higher energy 
range (1-10Bev). ‘These produced secondary particles in the atmosphere 
which were registered by neutron monitors as ground level effects (G.L.E.), 
bringing up to ten the number of such events which have been so recorded since 
1942 February 28. While the energies of flare accelerated protons are probably 
continuous over a wide range, it is useful to distinguish between the low-energy 
and high-energy events as P.C.A. protons and G.L.E. protons respectively. The 
flare accelerated particles are thus known to cover a wide range of energies 
(cf. Section 5). 

2. Class 3 flare of November 10.—The first of the Ha pictures of Plate 7, 
taken at og" 50™ 50° U.T. with the Lyot heliograph at the Royal Observatory, 
Cape of Good Hope, shows the appearance of the spot group 25 minutes before 
the flare began. The shape of the bright hydrogen plage to the east (left) and to 
the north of the two main umbrae should be noted, since it was in this region that 
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flare of 1960 November 12 showing the nimbus (dark halo), e.g. 14h 13m 353s. 


The Class 3+ 
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the principal flare activity occurred. At 10"15™46* brightening began in five 
separate small points of this plage and these quickly expanded and linked up to 
form the bright filaments A, B and C (Fig. 1). 
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Fic. 1.—Development of Class 3 flare emission filaments, 1960 November 10. 


As we can see from the light-curves given in Fig. 2 for the different regions, the 
increase in light intensity was most rapid between 10" 15™ and 10" 21™ which we 


regard as the flash phase of the flare. 


At 102 17™ 458 the bright dot D made its 
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first appearance over the main spot umbra and remained stationary there, with 
some fluctuations in intensity, up to the last observation at 12"31™. The four 
regions which were measured photometrically reached their first brightness 
maxima in the order B, A, D, C between 10" 19™ and 10% 24™., 
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Fic. 2.—Photometric light-curves in Ha for the four emission regions (A, B, C and D) of Fig. 1. 
(See Section 3 for definition of I/Ig.) The circles give the measured areas (uncorrected for fore- 
shortening) of the total emission in terms of the unit 1 square degree heliographic at the disk centre. 


Thereafter their intensities steadily declined until about 10%55™, when B, C and 
D again began to increase in brightness, reaching secondary maxima at various 
times between 11" 05™ and 112 30™., 

This second maximum of the flare also showed itself in the absorption of 
short radio waves traversing the D-region of the ionosphere. For example, the 
B.B.C. reception of Paris on a frequency of 6-175 Mc/s shows a double fade-out. 
The ionospheric absorption reached its first maximum at 10% 32™ followed by a 
partial recovery of signal strength, and then the absorption increased again with 
its greatest effect at 11"30™. All signals on frequencies up to about 9:5 Mc/e 
were affected. 

The solar radio noise records showed a notable steep rise of intensity on 200 
and 545 Mc/s starting at 10" 20™ and fairly great bursts on gooo and 3000 Mc/s 
starting at 10%12™, However, this event was surpassed in importance by a 
very great burst starting at 11" 16™. coinciding with the second maximum of the 
flare and the area increases of regions B and C. 

The filaments A and C showed no movement or drift relative to the sunspots, 
and the circular bright dot D, being located almost centrally above the principal 
umbra, was likewise stationary. The northern edge of B on the other hand, 
appeared to move northwards through 30000km in the interval 104%24™ to 
11"15™, mean velocity 1okm/sec. The measured (uncorrected) areas of the 
flare are given in Fig. 2 (open circles), the unit being 1 square degree heliographic 
at the centre of the disk. It can be seen that there were two area maxima, one at 
108 30™ and the other at 11"30™. The area at the second maximum was 14:5 
square degrees, and if we correct this for foreshortening by means of the sech 
relation*, h being 38° for the flare, we obtain a corrected area of 18-4 square 
degrees, or 890 millionths of the hemisphere. 


* A=Am sec h, where A.is the corrected area in square degrees, Am is the measured (i.e. 
projected) area in terms of the unit r square degree at the disk centre, and h/ is the heliocentric 
angle. The areas plotted alongside the light-curves are measured areas (Am). 1% square degree 
heliographic is equivalent to 48-5 millionths of the hemisphere. 
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Fic. 3.—Development of the Class 3+ flare emission filaments, 1960 November 12. 
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3. Class 3+ flare of November 12.—This flare began with a development 
in the plage at the east end (left) of the spot group very similar to that of 
November 10 but, unlike the former, the bright filaments extended into the west 
end of the group and covered up both of the large spot umbrae. 

The drawings of Fig. 3 and the photographs of Plate 8 illustrate the immensely 
rapid extension of the filaments A, B and C during the flash phase from 
13%23™ 358 to 13"28™ 358. In the interval 27™35° to 32™35° the filament B 
broadened to such an extent that it completely covered the two large spots and 
C flashed up suddenly at the same time. Much more frequent exposures than 
one per minute would have been needed to follow in detail these extremely rapid 
developments. In the 60-second interval between the exposures at 26™ 358 and 
27™ 35% filament B had extended westwards 60000km and eastwards 40000 km. 
If we interpret these extensions as movements of the luminous material we obtain 
minimum velocities of 1000 km/sec and 670 km/sec respectively. 
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Fic. 4.—Photometric light-curves in H« for the three emission regions (A, B and C) of Fig. 3. 
(See Section 3 for definition of I/I,). The circles give the measured areas (uncorrected for fore- 
shortening) of the total emission in terms of the unit 1 square degree heliographic at the disk centre. 


The light-curves of the filaments A, B and C (Fig. 4) show the rapid increase 
in Ha intensity during the flash phase, with B and C reaching their maxima at 
about 32™35° and A some 3 minutes earlier. B and C decreased rapidly in 
brightness with a slower decline setting in after 13"40™. This time also marked 
the change-over from the earlier condition of rapid chaotic developments to a 
state of quiescent fading. Our last photograph was taken at 14%31™ and the 
ending time given by the Sacramento Peak Observatory was 19"%22™, i.e. 
duration ~6 hours. The maximum measured area was 38 square degrees at 
13" 39™, or 42 square degrees (2040 millionths of the hemisphere) corrected for 
foreshortening (h=25°). 

The drawings of Fig. 3 (and Fig. 1) were made by projecting the heliograph 
films of the solar disk on to a white circle 11-5 inches in diameter* and drawing 
in the flare filaments in black. They provide a much clearer appreciation of the 
filamentary structure than is obtainable from the photographic prints in which 
the range of contrast is so much less. 

The light-curves given in Figs. 2 and 4 express the intensities in terms of the 
ratio I/I,, where J is the measured intensity of the filament and J, is the intensity 

* This scale is such that 1 degree heliographic at the centre of the disk is equivalent to 0-1 inch. 
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of the local Ha background, measured at a point well away from the flare and 
having the same heliocentric angle (4). ‘The methods of photometry and the 
standardization of the films are fully described in an earlier publication (7). 

While the measures of I/I, give a true indication of the light variations of the 
flare, as transmitted by the filter (band-pass o-7A), it should be noted that they do 
not enable us to deduce the central intensities of the Ha line for the flare in the 
absence of a line profile. Our previous spectrographic work indicates that the 
peak central intensity of the Ha line in a 3 + flare is of the order of 300 per cent of 
the continuum. 

The X-ray emission from this flare produced intense ionization in the D-region 
of the ionosphere. There was a short wave fade-out having a very sudden onset 
at 13"24™ (coinciding with the H« flash) with recovery to normal conditions at 
about 17". The B.B.C. Tatsfield Receiving Station reported that all short wave 
frequencies up to 28 Mc/s were affected by this absorption. At Dunsink a sudden 
enhancement of atmospherics (S.E.A.) was recorded on 24kc/s, beginning at 
13%25™, maximum at 13"33™ and ending approximately at 17". The 
H-magnetograph also registered a geomagnetic crochet with start at 13%24™ 
and maximum at 135 31™, 

4. The flare nimbus.—In a recent note (8) we drew attention to what is 
apparently a new phenomenon associated with some Class 3 and 3+ flares. 
This is a dark halo, some 300 ocookm in diameter, which begins to form 
round the flare region within a few minutes of the time of peak light intensity 
of the emission filaments. 

We have studied five examples of the phenomenon—the flare nimbus*, as we 
propose to call it. These occurred in association with the great flares of 1959 
July 16, 1960 April 1 (8), June 1 (g) and those of November 10 and 12 described 
above. For the first case our data have been obtained from a film taken with the 
Lyot heliograph of the McMath-Hulbert Observatory through the kind co- 
operation of Dr Helen Dodson-Prince. For examination of the last four we have 
used our own films taken with a similar instrument at the Royal Observatory, 
Cape of Good Hope. 

Although the appearance of the nimbus is less conspicuous on photographic 
prints than it is on the original 35 mm films, its general aspect can well be seen 
in the illustrations of Plates 7 and 8, e.g., Plate 7: 1960 November ro (115 43™ 158 
as compared with og" 50™ 508), 1959 July 16 and 1960 June 1; Plate 8: 1960 
November 12 (14"13™358 as compared with 13523™35%). In the drawings of 
Figs. 1 and 3, for the flares of November 10 and 12, the approximate positions of 
the outer boundary of the nimbus have been marked in (dashed lines) from the 
times of its first certain appearance. 

In Table I we summarize some characteristic properties of the nimbus, such 
as times of beginning, greatest visibility and maximum dimensions in relation to 
the flare times. The nimbus makes its first appearance at, or within a few minutes 
of, the time of maximum flare intensity, as shown by the light-curves. In one 
case only (June 1) was it seen before maximum. Its most conspicuous phase 
follows about 30 minutes after the flare maximum, and its dimensions then lie 
within the range 2—4 x 10°km or 4':5—9’. Estimates of the ending times are 
subject to much uncertainty, but we can say that the duration is > 1—2 hours. 

* O.E.D.: A cloud, halo or atmosphere, investing a person or thing. 
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TABLE I 
Date Flare Nimbus Radio emission Cmc-ray 
effect 
Type IV 
ae Mean of times 
jameter ‘5. six ead 
5 
1959 Class GLa. u.T. 10° km 200 Mc/s 
July 16 3+ B 2118 B 2134 B 2117 G.L.E. 
M 2132 M 2211 ~3°0 M 2159 P.C.A. 
E >2320 E >2300 E 2517 B 2250 
1960 
April 1 3 B_ 0840 B_ 0903 25 B 0847 P.C.A. 
M 0859 M 0937 38 M 0908 Bog45 
E 1320 E1100 E 1047 
June 1 3+ B_ 0823 B_ 0850 1°6 B_ 0837 PCA. 
M_  ogoo M og10 2°4 M 0859 (?) 
E1630 E1030 E 0937 
Nov. 10 3 B_ 1013 B_ 1021 2°1 B_ 1020(1116)* P.C.A. 
M roar M 1100 2'°9 M 1124 
E >1230 E >1230 E 1350 
Nov. 12 3+ B_ 1315 B_ 1332 2°3 B 1327 G.L.E. 
M 1332 M_ 1400 3°5 M — B1328 
E >1425 E >1425 E 1757 P.C.A. 
B 1354 


* Two distinct events corresponding to the double maximum of the flare (see Section 2). 
Notation: For flares: B, M and E signify the times of beginning, maximum brightness and 
end. 

For the nimbus: B, M and E signify the times of first certain appearance, maximum 
visibility and end of visibility. 

P.C.A.: polar cap absorption in the D-layer, produced by solar protons (100-400 Mev 
range). 

G.L.E.: ground level effect of solar protons (Bev range), recorded by neutron monitors. 
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Fic. 5.—Photometric intensity plots across the sunspot group in a N-S direction (see Fig. 3) 
before the appearance of the nimbus (dotted line) and at the time of its greatest visibility (full line). 
Diameter of original solar rmage=15 mm. 

5. Discussion.—Fig. 5 presents a photometric examination of the nimbus. 
Intensity curves were plotted for two film frames on November 12, measures 
being taken along a N — S line intersecting one of the main sunspots (see Fig. 3). 
The first frame (13 23™ 35*) records the conditions before the nimbus appeared 
and the second (13% 59™35°) relates to the time of its greatest visibility. ‘The 
vertical separation between the two curves represents a darkening due to the 
nimbus of the order of 5 per cent of the continuum near Ha, and the diameter is 
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at least 300000km. We can rule out the possibility which we first considered 
that the nimbus might be caused by non-uniform development of the emulsion 
(Eberhard effect) at points in close proximity to the brilliant flare filaments ; 
for the most conspicuous phase of the nimbus always occurs long after (~ 30 min) 
the filaments have passed their greatest brightness. The reality of the pheno- 
menon is undoubted. We are also struck by the fact that all five cases we have 
examined have been associated with two other comparatively rare flare pheno- 
mena: (1) the acceleration into space by the flares of streams of protons having 
energies in the cosmic-ray range, and (2) the emission by the flares of radio noise 
radiation having Type IV continuum. 

The continuum radiation of Type IV was discovered by Boischot, who also 
investigated its properties (10). It usually begins near the time of flare maximum 
in Ha, has a duration of several hours and is characterised by a smoothness of 
intensity quite unlike the other types of flare radio outbursts, or the noise storms. 
The maximum intensity is reached some 20-30 minutes after the start of the 
emission. ‘The size of this radio source on the Sun has been measured on a 
number of occasions with the Nangay interferometer on a frequency of 169 Mc/s 
and has dimensions of the order 7'—12’. There is evidence for an initial motion 
of ascent into the corona, but during the period of decay the source appears to 
be stationary and its dimensions practically constant. In order to account for 
such properties Boischot and Denisse have put forward the view (11) that this 
radiation is of the synchrotron type andoriginates in a cloudof relativisticelectrons 
moving in a magnetic field. The power radiated on frequency 169 Mc/s requires 
about 100 electrons per cm® having energies of the order of several Mev in a field 
of a few gauss. It is probable that these high-energy electrons and the cosmic- 
ray protons have a common origin in the flare discharge. 

It is clear that the characteristics we have found for the flare nimbus, namely 
its time of commencement, time of maximum visibility and dimensions (Table I) 
are closely comparable with those of the Type IV continuum radiation. We 
are bound to conclude that there is a direct relation between the two phenomena. 
In other words, we regard the nimbus as the optical counterpart of the postulated 
cloud of relativistic electrons. 
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Fic. 6.—Trapping of high energy particles in the coronal region above a flare. The particles 
are assumed to oscillate between points of reflexion located in the chromosphere. 


It is not obvious why these high-energy electons should have any visible 
effects at the chromospheric level. However, the electrons are not stationary: 
each is moving in a helical trajectory along lines of force and will perform oscil- 
latory, movements between reflexion points in the magnetic field. A typical 
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magnetic loop in the coronal region above a flare is shown in Fig. 6. Magnetic 
fields of this type are a characteristic feature of the coronal (5303) films, such as 
those taken at the Sacramento Peak Observatory. They likewise show up in 
the great loop prominences (Ha) which frequently form over an active region on 
the limb following a major flare. These regions may well act as magnetic traps 
for high-energy particles during their lifetime, which may be of the order of a 
few hours. Energy decay through collisions will be most rapid at the reflexion 
points in the chromosphere. By analogy with the outer van Allen radiation belt 
the flare and sunspot group may thus be surrounded by an “‘auroral zone’’ in 
which these collisions are most numerous. We do not, of course, expect any 
emission of light, but rather the reverse, an absorption of light caused by the 
excitation of neutral hydrogen with a corresponding increase in the number of 
hydrogen atoms in the second quantum level. 

For two of the flares we have been considering, 1960 April 1 and June 1, we 
have reported (8, 9) changes in the Ha striation patterns surrounding the sun- 
spot groups which took place near the times of maximum flare intensity. We 
attributed these changes to a rapid partial destruction of the chromospheric 
magnetic field brought about by the flash of the flare. Both these flares exhibited 
the dark halo. At present we are inclined to regard the two phenomena—the 
halo and the changes of striation pattern—as quite separate and distinct mani- 
festations of the flare process. The changes in the striation patterns are the more 
difficult to observe, requiring for their investigation continuous seeing of the 
highest quality combined with film exposures of suitable density. 
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THE DUST CONTENT OF THE SMALL MAGELLANIC CLOUD 
FROM COUNTS OF NEBULAE 


A. $. Wesselink 


(Communicated by the Radcliffe Observer) 


(Received 1961 April 13) 


Summary 


Counts of nebulae in and near the Small Magellanic Cloud show a clear 
deficit when compared with outside fields. Contrary to current opinion 
it is shown that a normal dust content is in accord with these findings. 





Introduction.—In 1951, H. Shapley discussed the question of transparency 
of both the Large and the Small Magellanic Cloud on the basis of counts of 
extragalactic nebulae. Though he finds very definite evidence of absorption 
in the Large Magellanic Cloud, the following statement summarizes his view 
on the Small Cloud absorption: ‘‘...an examination of the tabulated numbers 
and of Figure 2 shows that the Small Magellanic Cloud is essentially transparent ”’ 
(x). Further we read (2): ‘‘ An interstellar absorption in the inner section of 
the Cloud of two or three tenths of a magnitude is not out of the question’’. 
In a later paper (3) Shapley states: ‘‘ My evidence for the transparency of 
the Small Cloud should not be overvalued because of the conspicuous irregularities 
in some parts of the sky in the areal population of faint galaxies’’. However, 
we read (4) in a radio astronomical study of the Magellanic Clouds: ‘‘An 
important result from this survey is that the amount of gas in the two Clouds 
is about the same whereas the Large Cloud is believed to contain much more 
dust than the Small Cloud. Thus the ratio of gas to dust, which has often been 
held to be constant everywhere is very different in the two systems’’. Further 
we find (5): ‘‘...although there is evidence that the Large Magellanic Cloud 
population is more extreme (i.e. younger) than the Small Magellanic Cloud. 
The latter view is supported by the greater quantities of ionised hydrogen and 
dust associated with the Large Magellanic Cloud’’. Again (6): ‘‘...On the 
other hand, dust is observed in the Large Cloud but not in the Small Cloud’’. 

In spite of Shapley’s own doubt with regard to the transparency of the 
Small Magellanic Cloud as found by him in 1951, these quotations show clearly 
that the absence of dust is generally considered to be a well-established fact by 
others. 

In view of the importance of the question in relation to evolutionary problems 
of the Small Magellanic Cloud, we have tried to reassess the situation using new 
counts of nebulae. 

Observations.—The plates were taken at the Newtonian focus of the Radcliffe 
74-inch telescope. The scale is 1mm=22"-45. The emulsion was Kodak 103aO 








504 A. }. Wesselink Vol. 122 


almost without exception though in a few instances 103aE, 103aD have been 
used. As the faint extragalactic nebulae are red they were expected to show up 
more easily on the latter emulsions; moreover the extinction should be smaller 
at the longer wave-length. However, not much difference was found in the 
usefulness of the three emulsions for our purpose. It was rare that nebulae 
visible on the one brand could not be seen as well on the other kind if exposure 
times and quality were similar. The 74-inch telescope was used at its full aperture. 
The plates were of quarter-plate size covering 0:27 of a square degree. The 
standard exposure time was half an hour but exposure times ranging from 
15 minutes to one hour have been used. 

In an investigation of this nature it is all-important to discriminate between 
faint nebulae and faint stars. Accordingly the observations of a particular field 
were continued until a plate of acceptable quality was obtained. Though the 
larger part of the material used was specially obtained for the present problem, 
several plates already existing at the Observatory could be used to advantage 
and were searched for nebulae. These were especially found among the long 
series of plates taken of NGC 121, NGC 1466, the Sculptor system and a field 
00" 47.373 ° 46’ (1875), of a region in the SMC chosen for the study of variable 
stars. Clearly it was not difficult to find plates of excellent quality of regions 
where there were many plates. Apart from the material just mentioned 45 plates 
with a total exposure time of about 35 hours were taken for the present study. 
The somewhat random choice of fields within the main body of the Cloud and 
close to it is indicated in Table I. Of these, 28 plates with an exposure time of 
23 hours were considered of sufficient quality for the nebular search. 

The search for nebulae was usually made on pairs of plates using the blink 
microscope kindly loaned to the Radcliffe Observatory by the Royal Observatory 
at the Cape. 

In many cases only the best available plate of a region was examined but 
frequently two good plates were helpful for discrimination between plate defects 
and genuine nebulae. The plates were searched three times, whereafter a decision 
was taken on the number of genuine nebulae in the particular field together with 
an estimate of the quality of the result. The quality was given on a scale of 1-10, 
1 being very poor and 10 being perfect. The figure was determined almost 
entirely by the seeing. In spite of the fact that every effort was made to secure 
at least one plate for every field in good seeing, results with qualities ranging 
from 4-9 have been discussed. Quality 10 was never assigned. 

After completion of the counts, the results were rendered as homogeneous 
as possible by applying corrections dependent on quality, exposure time and 
atmospheric extinction. Variations in reflectivity of the mirrors, sensitivity of 
the emulsions, and the photographic processing were treated as accidental errors. 

Unfortunately our material was not extensive enough to derive the relation 
between count and quality independently. I have therefore used the corre- 
sponding relation as derived by Hubble in his great survey of faint extragalactic 
nebulae (7). Accordingly I have equated, admittedly a rough procedure, his 
FP with my quality 3, his E with my class 9. I further interpolated linearly 
for his other classes. At Pretoria maximum transparency does not usually 
coincide with best seeing. If, therefore other conditions prevail at Mt Wilson 
there is some further reason to query the accuracy of our procedure. However, 
we should realize that the results of a study like the present one are subject to a 
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variety of errors and that rough methods of reduction are both unavoidable and 
warranted. The counts were all reduced to quality 8. 

The reduction for exposure time was effected by means of the relation found 
empirically by Hubble (8): 


A log n= 1-29 A log (exposure time). 


The counts were reduced on an exposure time of 30 minutes using this formula. 

The formula for the correction for atmospheric extinction derives from the 
extinction coefficient for photographic light and the formula describing the 
variation of m with apparent magnitude in a uniform universe. For this reduction 
I have used Table VI from Hubble (9). The counts were reduced to the zenith. 

In Table I the results for the various fields have been collected. Columns 
1 and 2 give the position of the centre of the field for 1875. Column 3 gives 
the distance of the centre of the field from the position o® 50™-0, —73°15’, 
which is somewhat arbitrarily taken as the centre of the Small Magellanic Cloud. 
Column 4 gives the actually counted number of extragalactic nebulae 2,. 
Column 5 gives the number derived from column 4 but reduced to the standard 
conditions of 74-inch telescope, Kodak 103aO emulsion, 30 minutes exposure 
time, quality 8 and zenith. Column 6 gives the galactic latitude. 

It is seen that there are no fields with galactic latitudes smaller than 40°. 
As there seems to be very little, if any, galactic absorption in the direction of 
the Small Magellanic Cloud (10), we did not correct for it, though the 
difference between a field at the galactic pole and one at galactic latitude 45°, 
according to Oort’s discussion of the Hubble counts (11), would amount to 
o™-50—0™-31=0™19. The corresponding Alogz is 0-6 x o-19=0'12. 

In Fig. 1 the numbers n, have been plotted against the distances A. The 
dotted line indicates the mean n, for the eight outlying fields with A> 4°, which 
is n,=16-7. On the right-hand side of the figure a scale is shown indicating 
various amounts of absorption on the assumption of no Magellanic Cloud 
absorption to correspond with the above figure of 16-7. The other scale divisions 
for various absorptions A were computed from the relation 


I 
A= = (log 16-7 —log,). 


Accuracy of the results—The accuracy of the directly counted numbers n, 
is in the first instance controlled by statistical fluctuation which is measured by 
its standard error: +4+/(n,). The figure is further inaccurate as a consequence 
of misses which reduce the number. On the other hand the figure is erroneously 
increased by the counting of objects as nebulae which really are poorly resolved 
close stellar images. The inclusion of cloud nebulosity increases the count 
erroneously. Hubble is of the opinion that the misses and the positive errors 
will compensate each other on the average. 

The n, or the reduced numbers were derived from the m, by multiplying them 
with a factor depending on quality, exposure time and atmospheric extinction. 
The uncertainty of the n, therefore involves both that of the m, and that of the 
factor. 

Discussion.—In spite of the uncertainties affecting the n,, Fig. 1 shows a marked 
lack of nebulae in and near the Small Magellanic Cloud as compared with the 
numbers found outside it. If interpreted as due to absorption, our results seem 
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TABLE I 
1875 
a A Ne %, b 
h m °° ° e 
21 58-6 —67 07 15°3 18 9°3 — 44 
23 169 —72 19 4°7 25 13°8 —45 
00 OI'5 —72 54 #5 12 19°0 —45 
00 19°4 —73 29 2°3 14 17°8 —44 
00 2I°I —72 13 2°4 9 9°8 — 46 
co 338 —73 49 13 3 3°8 — 44 
co 39°8 —74 31 I'°5 2 3°2 —44 
©0 40°0 —72 §7 08 4 2°2 —45 
00 47°2 —72 27 08 2 I'l —45 
©0 47°53 —73 46 0°6 ° ro We) —45 
©0 §2°5 —74 24 12 7 3°6 — 44 
00 54°0 —34 22 38'9 59 29°0 — 83 
©0 54°9 —72 51 o'5 ° roe) —45 
00 58:0 —71I 31 1°9 2 3°9 —46 
00 59°4 —72 43 o"9 ° oo —45 
OI 00°4 —73 52 I'o I 06 —44 
OI 04°5 —73_33 I'l I I°3 —45 
OI 06°3 —73 38 2 3 2°1 —44 
o1 06°6 —73 25 $2 I o7 —45 
OI 10°0 —74 00 1°6 2 I‘o — 44 
OI 40°9 —44 06 29°7 14 28-2 —70 
02 03°6 —72 12 5°7 10 7°8 —45 
02 34°0 —35 02 39°5 22 12°3 —65 
03 45°2 —72 04 12°8 14 17°8 —40 
03 50°0 —68 50 14°8 II 15°5 —4I1 
20. 
t tT tT a 
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Fic. 1.—Numbers of extragalactic nebulae (ordinates) against angular distance of field from 
centre of Small Cloud (abscissae). 
The dotted line corresponds to the mean of 8 fields each at more than 4° from the adopted centre. 





to indicate an amount of rather more than a magnitude; an absorption of the 
order of o™-2 or o™-3 as advocated by Shapley (2) seems to be quite out of the 
question. Our results are therefore at variance with the belief that the Small 
Magellanic Cloud is transparent and there seems to be no reason to doubt 
that it might contain a normal amount of dust. 

Unfortunately it is not possible to give an accurate estimate of the absorption 
from our data. This is a consequence of the fact that the number of nebulae is 
not only reduced by absorption but also by the star density. In a transparent 
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but crowded field the number of extra galactic nebulae counted will be smaller 
than in regions with low star density. 

The examiner of the plates is less certain of the true character of a nebulous 
image in a crowded field, since he knows that poor resolution of close stellar 
images is frequently occurring there. The writer has carried out an experiment 
to prove this effect by superposing two negatives, film to film, and so counting the 
nebulae present on both plates. The one plate contained relatively few stars. 
but several genuine nebulae. The other was a negative of a dense field in the 
Small Magellanic Cloud. While the combination contained hardly more stars 
than the crowded negative, only about half the number of nebulae present on 
the other could be found. 

Our results depend on the ratio of the number of nebulae found in and outside 
the Small Magellanic Cloud. It is therefore differential and does not depend 
on the actual number of nebulae brighter than a certain limiting magnitude. 
It is for this reason that we did not find it necessary to correct for coma as was. 
done by Hubble in his great survey of nebulae. But Hubble was interested 
in the actual number of nebulae down to a certain limiting magnitude. It is 
nevertheless of some interest to compare our results, uncorrected for coma, with 
those of Hubble (uncorrected). 

Hubble (12) finds 92 nebulae per square degree with the Mt Wilson 100-inch 
telescope, 1 hour exposure, Eastman 40 emulsion, at the pole of the Milky Way, 
and at the zenith. 

In my experience Kodak 103a0 is about 1 magnitude more sensitive than 
Eastman 40. Reducing my result of 17 nebulae per 0-27 square degree, 
30 minutes’ exposure and Kodak 103aO emulsion and Radcliffe 74-inch telescope 
to Hubble’s conditions, I find 102 nebulae. Considering the many uncertainties 
involved in such a comparison the agreement between Hubble’s and my results. 
is surprisingly good. 


I am greatly indebted to Dr A. D. Thackeray who encouraged me to publish 
my results and who gave me the benefit of many discussions. Dr Thackeray, 
very generously, allowed me the use of many plates taken by him. 


Radcliffe Observatory, 
P.O. Box 373, 
Pretoria, 
South Africa: 
1961 April. 
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XUM 


ABSORPTION AND REDDENING IN THE MAGELLANIC CLOUDS 
A. }. Wesselink 
(Communicated by the Radcliffe Observer) 
(Received 1961 April 13) 


Summary 


The frequency distribution of colour excesses of stars brighter than a 
given apparent magnitude in distant stellar systems is studied. The fre- 
quency distribution (Fig. 1) depends on the luminosity function and the ratio 
of total absorption to colour excess. It is shown that most stars observed are 
only slightly reddened even in systems that are optically thick. There is 
no contradiction with the results of the preceding paper in which nebular 
counts near the SMC have proved considerable quantities of dust. 





1. Introduction.—In the preceding paper (1) we discussed the absorption 
in the Small Magellanic Cloud as shown by the reduction of the numbers of 
extragalactic nebulae. It was concluded there that total absorptions of over 
one magnitude may occur, at least in the more central regions, and that in these 
parts the system may be optically quite thick. 

In the present paper we discuss the information on absorption in a distant 
stellar system as obtainable from colour excesses. In this connection the frequency 
distribution of colour excesses of stars brighter than some apparent magnitude 
is of interest. . 

2. Frequency distribution of colour excesses of stars brighter than a certain limit 
of apparent magnitude.—Consider a certain area in a distant stellar system in 
which the optical thickness is more or less constant. Let the system be so distant 
that percentage differences in distance of the member stars from the observer 
can be neglected. 

It is then easily seen that provided we make a suitable adjustment of the 
zero point of apparent magnitude m we have: 

m=M+A (1) 
where M is the absolute magnitude and A is the absorption. 

Let the number of stars in the area having absorptions between A and 4+dA 
be g(A)dA; and suppose ¢(M)dM represents the number of stars having 
absolute magnitudes between M and M+dM. 4(M) is therefore proportional 
to the luminosity function. The number of stars in the area having M between 
M and M+dM, and A between A and A+4dA, provided there is no correlation 
between absorption and absolute magnitude, is then: 


$(M)q(A)dM dA. (2) 
The number of stars having m between m, m+dm and A between A, A+dA 
is then proportional to 


$(m—A).q(A)dmdA. (3) 
38 
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The brighter part of the luminosity function can generally be represented with 
fair accuracy by a form: 


¢(M)~e™. (4) 

Hence, adopting (4) we may write (3) 
e*™ .e-*40(A).dmdA (5) 
and the number of all stars having m between m and m+ dm is proportional to 
e*” dm (6) 


which follows from (5) after integration over A. 

The same exponential form as (6) holds for the number of stars brighter 
than some apparent magnitude and clearly is valid not only for a particular area 
but holds for the entire system. 

We now turn to the distribution function of absorptions A. At a given 
apparent magnitude the number of stars with absorptions between A, A+dA 
is proportional to 


e~*4.q(A) dA (7) 


and since this expression is independent of m it is valid not only at a particular m 
but for all stars in the area that are brighter than m. 

The luminosity function differs from one stellar system to another mainly 
as a consequence of differences in population. Even in one system variations occur 
from one region to another. It is therefore not possible to give a definite value 
for « which is valid generally. However, we are only concerned with the order 
of magnitude of « and shall accordingly adopt a plausible value. 

Schilt, Epstein and Hill (2) derive a value of Alog N/Am=o0-60 where N is 
the number of stars between m—4, m+} for a certain region of the Small 
Magellanic Cloud. This value 0-60 is close to van Rhijn’s result for the bright 
end of the luminosity function (3). de Vaucouleurs (4) finds Alog N/Am=0-76 
for the Small Cloud and only 0-43 for the Large Magellanic Cloud. We shall use 
Alog N/Am=o0-60 as a rough average value. The corresponding value for « is 
0°60 x log, 10= 1°38. 

If Amax is the maximum absorption in the area, in which the maximum 
amount of absorption is assumed not to vary much, q(A) is not likely to vary much 
in the interval o<A<Amax. This is exactly true for a system homogeneous in 
stars and dust. However, the actual absorption is more likely to be patchy. 
It is possible to calculate g(A) for such a system on the assumption that it is 
filled with randomly distributed clouds, each of which is capable of absorbing 
a small but definite amount. Again, the variation of q(A) with A is only slight. 
Since, then, there is no doubt that the exponential in (7) is the more varying 
and therefore more important factor we may take g(A) constant with sufficient 
accuracy and represent the distribution function of the A as being proportional to 


e184 dA (8) 


in the interval o< A < Amax with a frequency zero for values A larger than Amax. 
The function is shown graphically in Fig. 1. 

We might consider the average absorption for stars brighter than a given 
apparent magnitude using (8). Because the distribution function is so skew 
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stars selected according to apparent magnitude. 
7) a ; 
it is actually simpler to consider the median absorption. Denoting median 
m values with a ~ above the quantity we easily find in case of large Amax 
ily A= loge2 (9) 
ur S 
ue which, with «= 1-38, becomes 
a A=o"-50, (10) 
is It is of interest that this result is independent of the amount of absorption per 
all parsec. It depends on the luminosity function and the assumption of Amax 
ht being large. It follows that two systems which are both optically thick, which 
76 have the same luminosity function will have equal A even if their absorptions 
ise per parsec are different. However, the observer will see further in the system 
is with the smaller absorption per parsec. 
If the luminosity function used corresponds to photographic light and if the 
im selection of stars is according to photographic apparent magnitude, we derive 
ich a median colour excess on the B—V scale, E,_,=A/4:-0=0™-12 if we adopt 
in A, ,/Ez-y=4°0. When the stellar system is not optically thick we shall find 
ny. smaller values for the median colour excesses than the values just given. It 
is follows from the analysis that small amounts of reddening do not justify the 
ing conclusions that the systems must be transparent. 
ht. In a study of the brightest stars in the Magellanic Clouds, Feast, Thackeray 
ing and Wesselink (5) have discussed the colour excesses of highly luminous B stars 
ent in both Magellanic Clouds. They find an average reddening of about o™-10 
| to in either Cloud but it is not quite certain by how much this has to be reduced 
8 because of galactic reddening. It is, however, certain that some reddening occurs 
(8) amon in the Clouds th 1 Thi ll with 
g our stars in the ouds themselves. is agrees well with our 
ax- theoretical prediction for stars selected according to apparent magnitude even 
if the Magellanic Clouds were both optically thick. 
yen We therefore do not agree with Arp’s statement ‘‘that some SMC stars 
ew are unreddened and that the remainder are no more than slightly reddened ”’ (6). 
38* 
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Down to a certain apparent magnitude most stars may be only slightly 
reddened, but the few which have large colour excesses are of greatest interest 
for the problem of the optical thickness. The discovery of a single star highly 
reddened by dust within the Magellanic Cloud will give a lower limit to the 
maximal absorption in that direction. So far the greatest value for the visual 
absorption established in the Large Magellanic Cloud is o™-96 whereas the 
largest value found in the Small Cloud is o™-60. 

Though the discovery of moré stars in the Magellanic Clouds with larger 
reddening than the two just mentioned would be extremely valuable we may 
yet say that the present data on colour excesses are consistent with my results 
from nebular counts (1). It was concluded in that article that the Small 
Magellanic Cloud, like the Large Magellanic Cloud, contains a normal amount 
of dust. 


Radcliffe Observatory, 
P.O. Box 373, 
Pretoria, 
South Africa: 


1961 April. 
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” AND THE HEAT FLOW THROUGH THE LUNAR SURFACE 
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ger 
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ults 
nall Summary 
unt Observations of thermal radiation from the Moon at a frequency of 


178 Mc/s are described, giving a value for the mean disk temperature of 
233 +8 °K, which is the same to within 25 °K as various estimates of the 
mean disk temperature at the surface. The absence of any detectable steady 
temperature gradient through the surface leads to upper limits on the steady 
heat flow from the Moon for a range of possible properties of the surface 
materials. The observations at shorter radio wavelengths favour a model in 
which a gravelly material extends to depths of at least several metres. For this 
model the upper limit to the heat flow is very close to that expected from 
radioactive decay in a Moon of chondritic composition. 





1. Introduction.—Observations of the thermal radiation from the Moon 
in the infra-red and at radio wave-lengths have now been made by a large 
number of workers. The radio measurements show a variation of the apparent 
temperature with the lunar cycle which is much smaller than that observed in 
the infra-red and which lags in phase behind the cycle of illumination. A material 
having a very low thermal conductivity and a small absorption coefficient for 
radio waves has been proposed to account for these results since, for such a 
substance, most of the thermal radio emission originates at depths of several cm 
where the variation of temperature is much smaller than at the surface. The 
absorption coefficient must decrease with frequency since at lower frequencies 
no variation of temperature with the lunar cycle has been detected. 

A measurement of the apparent lunar temperature at a low frequency can 
therefore provide a value of the mean temperature far below the surface. If 
this differs from the mean surface temperature, then a knowledge of the probable 
value of the thermal conductivity of the materials forming the surface layers 
leads to an estimate of the steady heat flow through the lunar surface. 

In the present paper a measurement of the lunar temperature at a frequency 
of 178 Mc/s is described first of all. The results of previous observations at 
higher frequencies are then analysed to determine the probable values of the . 
thermal conductivity and penetration depth for radio waves of the material forming 
the lunar surface. Upper limits to the heat flow through the surface are derived 
for a number of models and are compared with the expected heat flow due to 
radioactivity in the Moon. 

2. Observations at 178 Mc/s.—The observations were made with the larger 
section of the interferometric radio telescope at Cambridge (1), which has‘beam 
widths to half-power points of 13’-6 in right ascension and 4°-6 in declination, 
used with a total power receiver. Sources are observed at meridian transit 
with this instrument and the procedure in the observations was to record a drift 
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curve across a region of sky containing the Moon on one day and a comparison 
record when the Moon was absent on the following day. Portions of two 
recordings of this kind are shown in Fig. 1. 
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Fic. 1.—Recordings of total power at 8=-+-07°. Full line: 13.11.60; dotted line: 14.11.60. 


At this frequency the sky temperature is comparable with that of the 
lunar surface. A determination of the lunar temperature, 7,,,, therefore depends 
upon an accurate knowledge of the sky temperature, 7’, in the position of the 
Moon but is relatively insensitive to the properties of the particular aerial system 
being used to measure the difference (7,,— 7 ,). ‘This feature of the observations 
contrasts with the microwave measurements where the chief uncertainty is in 
the forward gain of the aerial system. 

The results of several determinations of the lunar temperature at 178 Mc/s 
are assembled in Table I. 


TaBLe I 
Date U.T. of a ) AT, °K (T,,.—T,) <i i 
transit 
hm 3s w..¢ 
15.9.59 2304 22 40 44 —06 oo 57ot1'0 81+16 157 238+ 16 
16.9.59 2351 23 32 09 —OI 54 5‘OtI1'0 150+30 170 320 + 30 
17.8.60 0812 05 54 46 +18 20 OO+I'5 o+27 347 347427 
19.8.60 0947 07 37 34 +17 20 6-0+3°0 98 +49 180 278+49 
20.8.60 1034 08 28 29 +15 36 S$. +3 174465 107 281 +65 
13.11.60 0717 10 46 11 +07 31 8-4+0°8 118+10 117 235 +10 
27.11.60 10946 ©0 12 53 —00 34 IDo+1'5 203 +27 161 364427 


28.11.60 2034 OI 03 43 +03 46 6r+1'5 173 +42 186 359+42 


To convert the observed difference in aerial temperature with and without 
the Moon, AT, to the true temperature difference, (7,,,— 7,), account has been 
taken of three factors: 

(i) The fraction of the aerial beam occupied by the Moon. 
(ii) The position of the Moon in declination relative to the centre of the beam. 

(iii) The aerial efficiency or the fraction of the total power in the main beam of 
the aerial. Accurate measurements of the side-lobe pattern have given a value 
of 0-68. 

The quoted values of sky temperature, 7, were obtained from two series of 
observations made recently at Cambridge. 

(a) A survey of the region of sky —05°<5< +60° made by Turtle and 
Baldwin (in preparation) with the large 178 Mc/s aerial. 
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(5) Observations by Turtle (in preparation) at 178 Mc/s using an aerial 
of low gain similar to that used by Costain (2) in which particular attention was 
paid to obtaining an accurate zero level and sensitivity scale calibrated against 
thermal loads. 

The sky temperatures obtained in (a) were measured relative to the North 
Pole whose absolute temperature was not determined. ‘Turtle has smoothed 
the results of (a) to the beam width used in (b) at several different positions on 
the sky, thus both providing the detailed survey, (a), with an accurate zero level 
and also a check on the sensitivity scale. 

Some of the daytime records were affected by reception of solar emission in 
weak side-lobes of the main beam, in particular those of 1960 August 17, 19 
and 20. The remaining results show a scatter which is much larger than the 
probable error of each determination. Several effects might lead to such a 
variation. 

(i) Variation with the lunar cycle similar to that observed at wave-lengths 
near 1cm. Observations at 22cm by Davies and Jennison (3) have shown 
the variation to be <o-7 per cent and it should be negligible at a wave-length 
of 168 cm (v=178 Mc/s). 

(ii) The galactic radiation reflected from the Moon’s surface will vary with 
the position of the Moon in the ecliptic. The results of radar observations have 
been interpreted by Evans (4) as indicating a reflection coefficient at normal 
incidence of 0-06. For most positions of the Moon in the ecliptic the reflected 
radiation will contribute ~ 30 °K to the mean disk temperature but the variation 
in this quantity will be less than 4 °K except at times when the Moon crosses 
the galactic plane. 

(iii) Gradients of sky brightness across the aerial beam may lead to incorrect 
estimates of 7',, but in all regions of sky used the general gradients are very small. 
Occasionally, however, the Moon will occult sources of small angular size. There 
is a probability of 1 in 16 that a source of flux density 1-0 x 10-*6 w. m~*(c/s)-? 
is obscured by the Moon and in such an event the estimate of T7,,, is about 15 °K 
lower than the true value. 

(iv) Reflected radiation from the Sun and from terrestrial transmitters may 
be very important. The reflected radiation from the quiet Sun contributes 
<1°K to T,, but even comparatively weak activity on the Sun can produce 
measurable reflected radiation. 

A transmitter on the Earth radiating 100 kw isotropically over a hemisphere 
at a frequency in the 4 Mc/s bandwidth of the receiver would give an increase 
in the apparent lunar temperature of 60°K. Interference of this kind could 
originate from any place where the Moon is above the horizon. Since the 
frequency of 178 Mc/s lies in a band allocated to television transmission in the 
United States and in Europe, such transmitted powers may well occur at certain 
times of day. The most favourable time for lunar observations on the meridian 
at Cambridge is at about 08"00™ U.T. All American stations visible from the 
Moon have local times later than 02" 00™ when the majority will have closed 
down. Most European stations will not yet have begun transmissions. The low 
value of T,, obtained on 13.11.60 at o7%17™ U.T. bears out this expectation. 

Fortunately it is possible to distinguish recordings affected by reflected 
radiation from those on which only thermal radiation is present. In the latter 
case the Moon behaves as an effectively black disk showing only slight limb 
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darkening, whereas for reflection it behaves as a comparatively well-polished 
sphere, the greater part of the reflected radiation coming from an area only 
about 3’ arc in diameter at the centre of the Moon. The patterns drawn out in 
these two cases can be distinguished with the 13’ beam of the aerial, as illustrated 
by the recordings of 28.11.60 and 13.11.60 shown in Fig. 2. 


40°K 


S 





(a) 13.11.60 (b) 28.11.60 


Fic. 2.—Plots of the difference of aerial temperature, 4Tg, on successive days showing 
(a) thermal radiation from the Moon and (b) reflected radiation. The dotted curves are the 
theoretical ones for these two cases. 


Only two of the measurements (15.9.59 and 13.11.60) fit the pattern expected 
for pure thermal radiation. The two values are reasonably in agreement giving 
a weighted mean value for 7), of 236+8°K. To correct this to the true value 
of the mean disk temperature requires a knowledge of the dielectric properties 
of the Moon which are at present uncertain, but the reflection coefficient at normal 
incidence is probably close to 0-06. The mean sky temperature at 178 Mc/s 
is 260-270 °K, only a little higher than the apparent disk temperature, 7,,. 
It can be shown that, in this case, 7), is approximately 3 °K higher than the true 
disk temperature. We shall therefore adopt 233 °K as the mean disk temperature 
in the following discussion. 

3. Theory of the radio-frequency temperature.—Analyses by Wesselink (5) 
and Jaeger (6) of Pettit’s observations in the infra-red during a total eclipse of 
the Moon (7) led to the conclusion that the outermost layer of the surface consists 
of a material such as finely divided dust having a value of (kps)~!2= 1000 
(k=thermal conductivity, p=density and s=specific heat). The rate of fall of 
temperature during the umbral phase of the eclipse indicated that the underlying 
layers probably have a lower value of (kps)—"? of about 125, similar to that of 
dry gravel. The eclipse observations provide information concerning the 
uppermost few mm of the surface from which most of the heat radiated during 
the eclipse originates. The nature of the material at greater depths can at present 
be derived only from observations of the apparent microwave temperature 
through the lunar cycle. 

The problem of the apparent microwave temperature of a homogeneous 
solid whose surface experiences a sinusoidal variation of temperature has been 
treated in some detail by Piddington and Minnett (8). If the time variation of 
the surface temperature is 7,coswt then the temperature variation falls off 
with depth x as exp (— Bx), where B=(wps/2k)"*. The intensity of radio waves 
passing through the material falls off as exp(—ax) where «=470/cel? 
(o=electrical conductivity and e=dielectric constant). Then for such a 
substance the apparent microwave temperature is 


T; 
(1+25+ 
where 5= /« and tan¢d =8/(1 +8). 





ai cos (wt — ¢) 
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To calculate the expected variation of microwave temperature at the lunar 
equator Piddington and Minnett assumed that the surface temperature was 
determined during the lunar day by equilibrium between thermal radiation and 
the incident solar flux (i.e. T occos!4 A where A is the longitude from the sub-solar 
point) and was constant throughout the lunar night. More recently Jaeger (6) 
has made calculations of the variation of surface temperature on the lunar equator 
assuming only a constant flux of solar heat and radiation from the surface 
proportional to 7*. These more accurate calculations provide, in principle, 
a much more detailed test of various models of the lunar surface. At present, 
however, most observations have been made with beam widths comparable with or 
greater than the angular diameter of the Moon. The only two components of the 
temperature variation which are detected with any efficiency are the steady 
component 7, and that having the fundamental period of a lunar month, 
T, cos (wt — $)/(1 +25 +267)"?, since higher harmonics are approximately 
averaged to zero over the lunar disk and are also more rapidly attenuated with 
depth than the fundamental. 

Thus for most of the observations it is sufficient to Fourier analyse Jaeger’s 
curves of equatorial surface temperature to give 7, and 7,. The values 
corresponding to three of his models are given below together with those 
corresponding to Piddington and Minnett’s assumptions. We shall now compare 
these theoretical values with the observations of the variation through the lunar 
cycle, from which 6 may be derived, and with the observed values of 7>. 


\ 
Taste II 

Model TK 7K 7Ti/To 
Piddington and Minnett 249 134 0°54 
Jaeger (kps)—1/2 = 1000 210 160 0°76 

(kps)-1/2 = 1000 layer 0-45 cm thick on top 
of (kps)~1/? = 125 222 141 0°63 
(kps)-1/? = 125 ; 238 124 0°52 


4. Radio observations of the lunar temperature.—The measurements of the 
lunar temperature made so far are assembled in Table III. 

The numerical values of TJ, and 7,/7,(1+265+ 267)? in columns 3 and 4 
are those quoted by the observers. Some of the observations were made with 
beam widths slightly smaller than the lunar disk and give essentially the 
temperature of the central point of the disk. These are marked with a letter c 
in column 5. The remainder are thought to be mean disk temperatures. For 
the steady component 7, the theoretical difference between these quantities is 
small; Ty ean aisk = 2°97 mean central Since the temperature falls off with latitude 6 
only as cos". 

The effect on the fundamental variation of period 27/w of averaging over all 
latitudes is also small and may be neglected, but the averaging in longitude 
reduces the apparent variation by a factor of approximately 7/4 relative to that 
of a point on the surface. 

If the Moon is assumed to be composed of materials having a dielectric 
constant of 2:8 (corresponding to a reflection coefficient at normal incidence 
of 0-06), then, at wave-lengths where the sky temperature is effectively zero, 
the apparent central temperature is 0-94 of the true central temperature and 
the apparent disk temperature is approximately 0-88 of the true disk temperature. 
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Tas_e III 


Observed Corrected Corrected 
Observer A A T,/To(1 +28+28))% =, =, /To(1 + 25+ 28")? 
cm °K °K 

‘Tyler et al. (9) 0°75 160 o"10 188 o'13 
Salomonovich (10) 0-8 197 0°16 c 209 0°16 
Gibson (11) 0°86 185 o'19 c 196 o'19 
Piddington (8) 1°25 215 O°17 252 o'21 
Zelinskaya et al. 

(12) 1°63 224 0°16 ec a7 0°16 
Grebenkemper 

(13) 2°2 200 0°05 235 0065 
Troitskii et al. 

I 3°2 170 <0°07 200 <0'09 
Castelli et al. (15) 10 256 <0°025 © a <0°'025 
Akabane (16) 10 315 ° 370 
Mezger et al.(17) 21 250 <0-'02 294 <0°025 
Waak (18) 21 230 a 270 —_ 
Westerhout (19) 22 270 — ai7 —_ 
Davies et al.(3) 22 270 <{0°0035 c 286 {0°0035 
Castelli et al. (15) 23 254 — 298 — 
Denisse et al. quoted 

in (20) 33 208 — 244 _— 
Seeger et al.(20) 75 185 <o'1o 209 <0'13 
Baldwin 168 236 —_ 242 — 


These correction factors have been applied to the observations to reduce them 
all to central temperatures on the Moon corrected for emissivity. A crude division 
of the observations into those in which the central temperature is measured and 
those in which the mean disk temperature is the measured quantity has seemed 
adequate for the present discussion. The corrected values of T, and the ratio 
of the fundamental component of the variation to Ty are given in columns 6 and’7 
The ratio of the varying component to the steady component has been chosen 
for comparison with theory rather than the absolute value of the variation in an 
attempt to avoid errors due to uncertainties in the calibrations of many of the 
radio observations. 

5. Derivation of the penetration depth, 1/«, for different models.—We shall 

consider first of all two models in which the surface layers are homogeneous and 
have the following properties : 
ModelI: p=2gem™, s=o-2calg?°K-!, k=2:5 x10-%calcm™sec™ °K. 
Model IT: p=2gcem-*, s=o-:2calg-1°K-1, k=1-6x 10-*calcm~!sec— °K-. 
‘These correspond to Jaeger’s models with (kps)~/?= 1000 and 125. The values 
of 7,/T, for these models are given in Table II and using the relations in 
Section 3 it is possible to calculate 5 and « as functions of wave-length, A ; 
a is further related to the loss angle of the dielectric, A, by tan A=Ae!?a/2z. 
‘The computed values of 5, « and tanA are presented for Models I and II in 
Tables IV and V respectively. The dielectric constant has been taken to be 
2:8 (Evans (4)). The values of « and tan A for 0-15 cm wave-length were derived 
from Sinton’s (26) eclipse measurements. 

Values of tan A quoted by von Hippel (21) are given for a number of substances 
in Table VI. Unfortunately there seem to be no published figures for 
iron—magnesium silicates of which chondritic meteorites are largely composed 
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(Urey and Craig (22)) and which are supposed by several authors to form the 
Moon’s surface. It is of interest to compare von Hippel’s figures with those 
deduced from the observations. 


TaBLe IV 
Model I (kps)-!?= 1000 
Acm 8 a tan A 
O'15 — 46 190 107? 
o-75 3°6 O13 44.107 
08 2°8 0°16 5°8 107? 
0°86 2°3 0°20 7-7 107-8 
1°25 2°0 ’ 0°23 13°0 107? 
1°63 28 0°16 11°7 107* 
2°2 78 0°059 5°8 107? 
3°2 25°5 <0°084 <12 107? 
10 >21 <0°022 <9°9 107” 
21 221 <0'022 <20 107 
22 > 160 <0-'0029 <2°8 107? 
TaBLe V 
Model II (kps)-/?=125 
Acm 5 a tan A 
O15 — 5°8 240 107% 
0°75 2°9 2°0 10-7 6°7 107% 
o'8 2°2 26 107? 9°4 107° 
086 1°97 3°4. 107? 13° 2 
1°25 1°6 3°6 107? 20 107% 
1°63 22 2°6 107? 19 160° 
2°2 6°3 9°2 107? g'I 107% 
3°2 =4°3 <1°3 107? <19 107% 
10 217 <3°4 107% <15 1073 
21 S17 <3°4107% “at 30° 
22 >130 <4°5 10-4 <4°5 107% 
TaB_e VI 
Substance Wavelength in cm 
300 100 10 3 1'2 
Fused quartz 10-4 06 10-4 ro-* 2°5 1074 
Muscovite 2 107-* 3 10-4 
Sandy soil dry 10% 6 107% 4107* 
Clay soil dry 20 107? 1g10-* 1°3 107% 
Soda-silica glasses 107? 2-3 107? 
Al-Mg silicates : 3-30 1074 6-50 10-4 


The values of tanA for Model II in Table V are similar to those of common 
substances while those for Model I in Table IV are several times larger than the 
highest values of tanA in Table VI. Any model assuming lower values of the 
thermal conductivity than 2-5 x 10~*calcm~!sec-! °K-! would require even 
larger values of tanA which would be scarcely acceptable. 

Unlike the values of tan A in Table VI, which show only a slow variation with 
wave-length, the computed values for Models I and II vary very rapidly with 
wave-length at wave-lengths in the neighbourhood of 1 cm. If this feature is 
regarded as real, then the electrical conductivity actually increases with 
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wave-length in the range A<1-25 cm. At wave-lengths greater than 1-25 cm, 
tan A shows a steady decrease with increasing wave-length. 

The calculation of the penetration depth, 1/x, at 168cm wave-length 
depends on the behaviour of tanA at longer wave-lengths than are included 
in Tables IV and V, but a safe lower limit to 1/« is obtained by assuming 
tanA to remain constant in the range 22cm<A<168cm. The penetration 
depths for the two models are then: 


Model I (1/«),16g em > 26 metres, 
Model II = (1/&) 45 om > 168 metres. 


Models in which both the thermal and electrical conductivities vary with depth 
present a very large number of possibilities and we shall consider only one type 
of model which seems physically the most probable. 
Model III; the composition of the surface is supposed constant with depth 
but the packing of the material increases downwards due to the pressure of the 
overlying material. Thus k may increase from about 2:5 x 10~* cal cm~! sec! °K-4 
at the surface to as much as 4x 10~%calcm~'sec~! °K~-! when solid rock is 
reached. The electrical conductivity on the other hand would be expected to 
increase by a factor only of the order of 2, depending on the relative density 
of the loose surface material and the solid rock. Using the relations for a 
homogeneous model from Section 3, 
— yj2 (PE)? a0 
§=B/a= 5 (ws/2)! a —— 

For a homogeneous model, 5ccA if tanA is constant. For the type of model 
considered here, the factor p"e!/?/tan A will vary but little with depth. At the 
shorter wave-lengths the low value of k at the surface is the relevant one, whilst 
at the longer wave-lengths which penetrate more deeply the effective value 
of k is much larger. Thus the apparent value of 5 should increase with 
wave-length more slowly than A!®. In fact, for A>1-25 cm, the observations 
show an increase of 6 faster than A'®. To account for the observations 
satisfactorily, tanA must decrease with increasing wave-length and there 
cannot be any significant increase of k with depth in the range (1/«),.o5em to 
(1/)o0 em: Below this latter depth we have no information regarding the value 
of k 

6. The steady temperature T,.—Piddington and Minnett’s and Jaeger’s 
treatments of the lunar temperature considered the solar flux as the only source 
of heat and thus the mean temperature, 7), should be the same at all depths in 
the lunar surface. A steady outward flow of heat, Q, due to radioactivity in 
the Moon, or a possible initial high temperature at the time of the Moon’s 
formation would give rise to a temperature gradient O/k °K cm~ at the surface 
and the apparent lunar temperature at a wave-length where the penetration 
depth is 1/x would be 





) ae ™~ 
T+ | Ox =e ax) dx 
0 


If « and k are independent of depth this becomes 


To + O/ka. 


We may now examine the observations to find any possible variation of JT) with 
wave-length and to check the values against the theoretical ones of Jaeger. 
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The spread of results in Table III reflects the poor accuracy of the calibrations 
of many of the observations. The value of T, at 168 cm lies about 10 °K below 
the mean of all the other observations, and thus there is at present no observational 
evidence to suggest that the mean temperature at 168 cm is systematically 
higher than at short wave-lengths. 

The values of T, in Jaeger’s models quoted in Table II range from 210-238 °K. 
Lower values of Ty could be expected only by assuming the surface layers to have 
a thermal conductivity <2:5 x 10-*calcm-!sec"!°K-!. This is lower than 
that measured for dust in vacuo and would also lead to improbably large values 
of tanA to explain the centrimetric observations. Values of 7) higher than 
238 °K would require a comparatively large value of k isemiilietail below the 
surface layer which Jaeger showed to be incompatible with the eclipse observations 
in the infra-red (Pettit (7)). He found that the model fitting Pettit’s results best 
was one in which a thin layer of dust lay on top of a gravelly base (corresponding 
to the third model in Table II). The mean temperature of this model is 222 °K. 
We may therefore conclude that the observed value of 7, at 168 cm does not 
differ from the mean surface temperature by more than 25 °K. Then 


O/ka<25 °K. 

For Model I 

O<2:5 x 10-*calcm~*sec™}, 
and for Model II 

O<2°5 x 10~* calcm~*sec—, 
Most models of type III, in which a layer of dust overlies a gravelly base, give 
intermediate values for the upper limit to Q. If, however, the thermal 
conductivity increases suddenly at depths just below (1/«)s. em to a large value 
(such as would be the case if solid rock were reached at this depth) the value 
of 7, at 168 cm would be approximately T.,.)+ O/Rag:em and the upper limits 
to O would be relaxed to give O<1-7x 10~*calcm~*sec~! for fine dust lying 
on top of the rock, and O<17x 10~*calcm~*sec™! for gravel lying on top of 
the solid rock. 

The steady heat flow through the surface for a Moon of chondritic composition 
has been estimated by Jaeger (23) as 2:3x10~‘calcm~*sec-! and more 
detailed calculations by MacDonald (24) gave values of 2-4 x 10-7 calcm~*sec~! 
for an initially cold Moon formed 4:5 x 10° years ago and 2-9 x 10~‘ calcem~* sec! 
for a Moon having an initial temperature of 600 °C. 

The upper limit on QO for Model I is only a tenth of the expected heat flow. 
This fact, together with the very large derived values of tan A, make this model 
unacceptable. A model of the lunar surface of the type proposed by Gold (25) 
in which large areas are deeply covered in dust can only be made to fit the 
observations if the pressure some distance below the surface increases the thermal 
conductivity by at least one order of magnitude. It seems doubtful if such a 
picture would then be consistent with the apparent variation of temperature 
through the lunar cycle discussed in Section 5. 

7. Conclusions.—(i) An accurate measurement of the mean disk temperature 
of the Moon of 233° + 8° K has been obtained at the lowest frequency which has 
so far been used in such measurements. 

(ii) A model of the surface consisting of fine dust of very low thermal 
conductivity extending to great depths seems unsatisfactory. 
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(iii) The absence, at 22cm wave-length, of any detectable variation of 
temperature through the lunar cycle indicates that gravelly material extends to 
depths of at least (1/«)22em=20 metres. 

(iv) For models which fit the centimetric observations best the observation 
at 168 cm provides an upper limit to the heat flow through the surface which 
is very close to that expected for a Moon of chondritic composition. 

To extend this investigation three types of observation are now needed. 

(i) Eclipse observations at the shorter wave-lengths to determine more 
precisely the appropriate model of the outermost layers which largely determine 
the mean surface temperature. 

(ii) Observations in the range 3-21 cm designed to detect any variation 
with the lunar cycle, preferably by comparison of the brightness temperatures. 
of the two halves of the lunar disk. 

(iii) An investigation of the dielectric properties of stony meteorites in the 
radio-frequency range. 
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discussion of the observations, and A. J. Turtle who has computed the accurate 
values of sky temperature used in the derivation of the lunar temperature. 
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ON THE CENTRE-LIMB VARIATION OF GRANULE CONTRAST 
R. G. Giovanelli 
(Received 1961 April 17) 


Summary 


It is shown that the observed centre-limb changes in granule contrast, and 
particularly the disappearance of granulation near the limb, are due mainly to 
foreshortening and finite telescope resolution. Without greatly improved 
resolution, it is invalid to use such observations for inferring the optical depth 
of the top of the granulation zone, as has sometimes been done. 





1. Introduction.—A number of recent papers (Résch 1957, Edmonds 1960, 
Loughhead and Bray 1960) have discussed the appearance of granulation near 
the Sun’s limb. The significance of the results has been in their implication as 
to the structure of the convection zone, at the top of which the granules apparently 
lie. Using a model in which the granules are supposed covered by a thin layer 
of optical depth 7, which reduces the granule contrast progressively from centre 
to limb, Plaskett (1955) related 7 to 0, the angle from the centre of the Sun at 
which granulation disappears (i.e. the contrast falls below some arbitrary level). 
Waldmeier (1939) and ten Bruggencate (1939) observed @ to be 70° and 72° 
respectively, corresponding.to which Plaskett inferred r= 1-2. Edmonds (1960) 
also found a rather similar result, about 75°, though he points out that a larger-scale 
lower-contrast pattern of brightness fluctuations seems to extend to 0=86° or 
more. On the other hand, Résch (1957) considered the granules to be visible 
to about 82°-85°, a result with which Loughhead and Bray (1960) are in 
agreement; de Jager (1959) puts the corresponding optical depth r=o°1. 

In the various discussions of these observations, no account has been taken 
of the contrast transfer functions (frequency responses) of the telescopes used. 
While the contrasts of very coarse structures are faithfully reproduced, they are 
reduced for finer structures, reaching zero at the resolution limit. Hence 
foreshortening must cause substantial reduction in the observed contrast towards 
the limb. The question arises as to whether this is a significant factor in the 
observations. As will be shown, instrumental effects have almost entirely 
dictated the results in all observations made to date. 

2. The contrast transfer function.—For simplicity I shall use the contrast 
transfer function for a one-dimensional structure; this describes the contrast 
with which the optical system reproduces each of the various Fourier components. 
in the object. Though granulation has a variation in two dimensions across the 
surface, it is sufficient for the present discussion to consider radial variations 
alone and to represent the granulation by a sinusoidal structure whose period 
corresponds to the granule diameter. 
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The contrast transfer function of an aberration-free objective may be written 
(e.g. Steel 1953a) 


T(m)= = [cos—!m — m(1 — m?)4?}, (m<1) ae 


) , (m>t1) 


where m=A/D4, A being the wave-length of light, D the diameter of the objective 
and ¢ the angular size corresponding to the spatial period. The function T(m) 
has been tabulated by Steel (1953b). A one-dimensional sinusoidal structure 
(parallel to the limb) which would have angular size ¢, at the centre of the disk 
is foreshortened at @ to 


$= $y cos 8. 


Fig. 1 shows the relative contrast of this structure (the relevant value of 7T(m)) 
as a function of 6 and of m.-1=4,)D/A. As guide, the values of ¢, corresponding 
to those of m,~1 in Fig. 1 are given in Table I for A=5-5 x10-5 cm and for 
diameters of the telescopes used by some of the authors. 


TABLE I 


Relationship between mo (Fig. 1) and do 
aad 


Mo 
1°5 2 3 5 10 15 20 30 
D (cm) Corresponding ¢» (sec of arc) 
12°5 1°36 18 2°7 4°5 gl 13°6 18-2 27°3 
23 0°74 1°0 1°5 2°5 570 74 9°9 14°8 
30 0°57 90°76 1°14 19 3°7 5°7 76 11°4 


Real telescopes are never aberration-free, while the solar image is also affected 
by atmospheric seeing. Aberrations do not influence the values of 7T(m) near 
m=o (very coarse patterns) or very near the limit of resolution, but cause a 
general lowering of 7T(m) for intermediate frequencies, or even a reversal of 
contrast for some frequencies in extreme cases; Blackwell, Dewhurst and Dollfus 
(1959) give the measured contrast transfer function of one telescope. The biggest 
reduction in T(m) is usually at about half the resolution limit. The effect of 
aberrations would be to lower all the curves in the central range of contrasts in 
Fig. 1. , 

3. Discussion.—Lacking knowledge of the actual frequency respunses of the 
various telescopes, I shall use Fig. 1 for discussion purposes. 

The diameters of photospheric granules generally lie in the range 1”-2”. 
For the largest of these (¢)=2”), observed with a 12-5 cm aperture telescope 
such as that of Loughhead and Bray, m,-1=2-2. From Fig. 1, the contrast has 
dropped almost to zero by 6=60°, beyond which such granules are below the 
limit of resolution. For Résch’s and Edmonds’ telescopes (apertures 23 cm 
and 30 cm) the corresponding values of @ are about 72° and 75° respectively. 
These conclusions are very little influenced by the use of a one-dimensional 
rather than a two-dimensional anaylsis. Edmonds’ observed values agree 
remarkably well with expectation, but the structures observed by Résch and by 
Loughhead and Bray extend far too close to the limb to be ordinary granules. 
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Fic. 1.—Variation in contrast transfer factor T due to foreshortening of a sinusoidal object at 
angular distance 8 from the centre of the disk, a single element subtending angle $y when at the centre 
of the disk. The curves are for various values of my-1=¢ D/A, D being the telescope diameter. 
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(a) Photospheric granulation photographed by Schwarzschild: telescope aperture 30 cm. 
(6) Photospheric granulation photographed by Loughhead and Bray: telescope aperture 12°5 cm. 
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(a) | (b) 


Photographs of Plate 9 (a) simulating observations of granulation through telescopes of apertures 
23 cm (a) and 12:5 cm (b). The direction of view varies from the normal as marked, showing the 
variation in appearance of granulation from centre (top) to near the limb (bottom) due solely to 
finite telescope resolution. The photographic processing of all the prints within either set (a) or set (b) 
is identical, but set (b) has been printed on paper more contrast than set (a). 
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It is possible to understand the various observations better by making use of a 
very simple demonstration. Plate g(a) is a photograph of the photosphere 
obtained at a high altitude by Schwarzschild with a 30 cm telescope. This has 
been rephotographed at various angles to the normal, at such distances and 
apertures as to simulate the centre-limb variation in granulation as viewed 
through a 23 cm aperture telescope (Plate 10(a)) and one of 12:5 cm aperture 
(Plate 10(5)). The contrasts of the finer structures must be somewhat lowered in 
the course of this processing. However, the upper photograph of Plate 10 (5) 
closely resembles the direct photograph of the centre of the disk, Plate 9 (5) 
obtained by Loughhead and Bray with a 12-5 cm aperture telescope. 

The general superiority of the larger over the smaller aperture is at once 
apparent, as is the rapid decrease in contrast towards the “ limb ” in both cases. 
The numbers of granules visible also decrease rapidly towards the “ limb ”’, 
groups of granules well defined at the centre of the disk merging into one apparent 
granule when foreshortened. A good example can be seen by noting the bright 
granules at the intersection of the arrows in Plate 10. Just where granulation 
“disappears” will obviously depend greatly on the criterion used by the 
individual observer; as pointed out by Edmonds, larger-scale brightness 
fluctuations, or rather groups of brighter granules, can still be detected when the 
individual granules are no longer resolved. 

It is clearly inadmissible to use the disappearance of granulation towards the 
limb for an analysis of the optical thickness of any layer above the convection 
zone, without much higher resolving powers than have been available hitherto ; 
or alternatively without a very much more sophisticated analysis of the variation 
in contrast with 6. An adequate analysis seems in any case impracticable without 
much higher resolution. The most that can be said at present is that any such 
layer cannot have an optical thickness exceeding about 1-2; it may very well be 
negligible. 


Division of Physics, 
Commonwealth Scientific and Industrial Research Organization, 
Sydney: 
1961 April. 
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OF GALACTIC RADIO EMISSION FROM A 
SATELLITE-BORNE RADIO TELESCOPE 
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Summary 


It has already been proposed to extend measurements of the spectrum of 
radio emission from the galaxy to frequencies of about 1 Mc/s by means of a 
satellite-borne receiver. If the satellite is in the upper ionosphere, where 
the refractive index is increasing with height, a focusing effect will occur 
which may allow individual features of emission to be observed. It is shown 
that a beamwidth of the order of 20° may be achieved under favourable 
circumstances; the most important limitation is the frequency bandwidth 
necessary in a sufficiently sensitive receiver, which inevitably means that 
dispersion in the refractive index effectively blurs out any narrower beam or 
any interference effects inside the beam. 

The experimental conditions needed for the veitoeion of this degree of 
focusing in a satellite or a rocket flight are discussed. It appears that the 
frequency of a satellite-borne receiver should sweep over a range of about 3 
to 5 Mc/s, with a bandwidth of 1okc/s. The satellite orbit should be 
nearly circular, at a height of about 400-500 km. 





1. Introduction.—Although theories of the galactic radio emission which 
attribute it to synchrotron radiation from relativistic electrons (Shklovsky 1952, 
Tunmer 1959, Hoyle 1960) show that the radio spectrum is not far from that 
derived from the cosmic-ray energy spectrum, there remains a discrepancy 
which is hard to explain. Briefly, it is expected that in the frequency range 
20-400 Mc/s over which accurate measurements have been made the radio 
spectral index should be near — 1-0, whereas the best radio measurements give 
an index of —0-37+0°04 (Costain 1960). A low index would only in fact be 
expected at lower frequencies, where the electrons lose energy primarily by 
ionization rather than by the synchrotron emission itself. Observations do not 
so far give any evidence of a change in spectral index with frequency, although 
it seems likely that an extension of the frequency range will show such a change. 
With observations available over an extended range of frequency it may become 
possible to revise the theory, and learn more about the processes of energy gain 
and loss in the cosmic-ray electrons. 

It is already proposed to extend observations below the frequency limit at 
present set by the ionosphere, by mounting a receiver for frequencies around 
1 Mc/s in a Scout satellite. A dipole aerial of about 150 ft overall will be used, 
with a receiver of about 20 kc/s bandwidth sweeping over a frequency range 
of about # to 3 Mc/s. No angular resolution is required for the purposes of this 
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experiment, which is intended only to provide a measure of the average brightness 
of one hemisphere of sky at any one time. However, there are several individual 
features of the galactic emission which deserve study with a finer angular 
resolution, and a beamwidth of the order of 10°-30° would be very valuable 
for this. With such a beam it would be possible for example to delineate the 
extent of absorption by Hu near the galactic plane, or to obtain a rough 
comparison of the spectra of such features as the spur of emission at /7=o, 

Such a beamwidth is impossible to obtain by ordinary means when the 
wave-length is of the order of 300 metres, but it will be shown that the focusing 
effect found in the upper ionosphere may provide such a beamwidth. It is likely, 
in fact, that this focusing effect will be observed to some extent in the Scout 
experiment, since the technique of frequency-sweeping is the most likely way 
of demonstrating it. 

This paper analyses the conditions under which focusing may be observed 
from satellites and from rockets. Observations of focusing are restricted by the 
complications of magneto-ionic splitting, ionospheric irregularities, and the 
extremely dispersive nature of the refractive index, but the practical examples 
show that a beamwidth of about 20° may be obtained. 

2. The focusing effect.—Ignoring for the moment the effects of the terrestrial 
magnetic field, we consider the upper ionosphere as a smooth horizontally 
stratified medium in which the refractive index increases smoothly from zero to 
unity in a height of some hundreds of kilometres. Rays to a receiver at a height 
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Fic. 1.—Ray path from satellite in the ionosphere. 


where .=/, received with an intensity J(«) in a direction «, may be traced 
back through the ionosphere until they emerge at an angle ¢ to the vertical 
(Fig. 1), where: 

4 sin d = py SiN a. 


The maximum value of ¢ is given by sin¢dy=py. All rays are contained in a 
cone of semi-angle ¢, so that a restricted aerial beam is obtained like that from 
mounting a dipole over a reflecting surface. The intensity J(¢) in direction ¢ 
is related to J(«) by: 


I(¢) sin ¢ dé = I(«) sin « dx 
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so that we may obtain the gain over an isotropic receiving aerial in direction ¢: 


cos p 


in?d\l2- 
(: _ sin *) 
Ko 


Here we have put J(«)=1, although practical aerials will depart somewhat 
from an isotropic sensitivity. The forward gain, at =o, is seen to be 1/9? 
and the gain increases towards infinity at $6=¢,, where there is a cusp of rays 
corresponding to values of « near go°. Fig. 2 shows a plot of the polar diagram 
I(¢) for py = 0-2. 
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Fic. 2.—Polar diagram I(¢) deduced from geometrical refraction. 


This simple picture of refraction suggests that narrow beamwidths can be 
obtained with a loop aerial system no more elaborate than a single dipole or loop. 
Subsequent analysis must take account of interference effects in this refraction, 
and of the effects of the Earth’s magnetic field. The ionosphere may be irregular, 
as is suggested by observations of the scintillation of radio stars; this may set 
a limit on the useful gain to be obtained. The choice of frequency must be 
considered carefully, since the refractive index changes very rapidly with frequency 
when it approaches zero. Finally, it may be difficult to attain the necessary 
sensitivity in the receiver, so that the choice of aerial system must be considered. 


a 








Fic. 3.—The two ray paths from the satellite. 
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3. Interference effects.—Using the simple model of the ionosphere it appears 
that for each angle of emergence ¢ there are two rays of equal intensity at the 
satellite, symmetrically disposed about the horizontal (Fig. 3). 

The path difference evidently depends on the form of variation of » with 
height, but the general effect is to provide a system of Airy fringes across the 
polar diagram. At the centre the order of the fringes may be quite low. For an 
ionosphere whose electron density varies as e~”*, it may be shown that the 


order m, at the centre is 
sry 2 
n= Hote eee] 


% L3 5 


where jy is the refractive index at the satellite, and Ay is the free space wave-length. 
For example, taking 4) = 100 km, Ay = 300 m we find values for mp) as follows: 


Ho orl o-2 03 O°4 
Ny 0°45 3°5 12 120. 


However, since ply is critically dependent on frequency, the order of interference 
will change so rapidly with frequency that it will be difficult to observe interference 
with a reasonable receiver bandwidth. For a small change in frequency Af 
we have a change in mp given by: 


Atty _ 3Apo _ , 1—Mo* Af 
— = == 3 z = 
No Ho bo” f 


and, for example, if 4»=0-2 a $ per cent in frequency gives a change An) ~1. 

Interference effects of this nature will not be considered further in this paper. 
A more detailed analysis of similar effects has been made by Haselgrove, Haselgrove 
and Jennison (1961), and by Budden (1961). 

4. The effect of the Earth’s magnetic field.—In the upper F-region the collisional 
frequency is low, and absorption effects may be neglected. The main complication 
in comparing the simple model with the actual ionosphere is therefore the Earth’s 
magnetic field. The refractive index is now different for the two characteristic 
modes of propagation, and for different directions of propagation. Fig. 4 shows 
the conditions under which refractive index is zero for propagation in directions 
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Fic. 4.—Magneto-ionic splitting of refractive index. The graph shows the frequency where 
#=0 at various electron densities (expressed as plasma frequencies f ,) for longitudinal and transverse 
propagation. 
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longitudinal and transverse to the magnetic field, for a gyrofrequency of 1 Mc/s 
(H =0-36 gauss). 

It will be seen that, at a fixed electron density, for all frequencies above a 
value f, on the line AB the refractive index is greater than zero, while at this 
frequency f, it is zero for all directions of propagation. This refers to the 
extraordinary mode of propagation; the refractive index for the ordinary mode 
remains positive at the frequency f,. The change in the refractive index for the 
extraordinary mode for frequencies slightly above f, is illustrated in Fig. 5 for 
the conditions f,=1 Mc/s, f,=1 Mc/s (Fig. 5(a)) and for f,=1 Me/s, 
fp=3 Me/s (Fig. 5 (6)). 
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Fic. 5.—Refractive index near zero. (a) fq=1 Mc/s, fy=1 Mc/s; (6) fa=1 Me/s, f,=3 Me/s. 
Upper curves refer to transverse, lower to longitudinal, propagation. 


It will be noted that the change in » with frequency is very rapid near the 
critical frequency, and also that above this frequency p» is dependent on the 
direction of propagation. If for example it were required to place a receiver 
where » =0-2 in the situation of Fig. 5 (a), the frequency could be adjusted so 
that w lay in the range o-1r§ to 0-25 according to direction (f= 1-63 Mc/s), but 
a receiver bandwidth of 10 kc/s would further widen this range to, say, 0-15 
to 0-30. At the higher frequencies of Fig. 5 (5) the situation is easier, with 
splitting between 0-20 and 0-26, widening to between o-20 and o-28 for a 
bandwidth of 10 kc/s. 

Splitting of ~ according to azimuth affects the shape of the beam in a way 
which depends on the inclination of the magnetic field. An example has been 
calculated by N. Capon, of the Cambridge University Mathematical Laboratory, 
for conditions where the inclination is high, giving only a small distortion of the 
beam. The model used for the ionosphere had a penetration frequency of 
6 Mc/s, and the observing frequency was 3-6 Mc/s. To obtain a refractive index 
Jp =0°2, the receiver was considered to be 270 km above the maximum, where 
the electron density was decreasing roughly exponentially with a scale height 
of 70km. The magnetic field was taken as a dipole field, and allowance was 
made for the curvature of the Earth. Fig. 6 shows the results of ray tracing for 
a position over southern England (magnetic latitude 55°). It is seen that the 
distortion of the beam is not serious. 

5. Irregularities—Local variations of refractive index are known to exist 
in the upper F-region. The condition of ionospheric echoes known as 
‘*Spread-F’’ indicates that the variations of density may often exceed 
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Fic. 6.—Angles of emergence ¢ for rays leaving a satellite at various elevations « for azimuthal 
directions. Observing frequency 3°6 Mc/s, scale height 70 km. Magnetic latitude 55°. 


one-tenth of the mean density. However, this condition is generally only 
observed at night time, and is even then rare in geomagnetic latitudes between 
about 15° and 45°. Clear F-region echoes correspond to less than 1 per cent 
variation of density; we may reasonably expect that less than 2 per cent 
variations usually occur. For u»=0-2, this corresponds to a variation of about 
0°05 in the value of pio. 

6. Practical limits of focusing.—The detection of galactic radio noise at these 
low frequencies with simple receivers suitable for mounting in satellites demands 
the use of appreciable bandwidths. Both from convenience in construction and 
from the necessary smoothing of the noise output when a short output time 
constant must be employed, it seems reasonable to take 10 kc/s as a typical 
bandwidth. From Section 4 we see that at 1 Mc/s this places a lower limit at 
fo =0°2 for useful focusing, where the frequency bandwidth together with the 
magneto-ionic splitting provides a maximum range of py) between 0-15 and 0-30 
according to direction. At 3 Mc/s the same bandwidth may be used to provide 
a maximum range of 4) between 0-09 and 0-17 or between 0-20 and 0:28. Clearly 
focusing may be achieved only down to yy~ 0-15 or 0-2, the changes being better 
at higher frequencies, and at high geomagnetic latitudes. This limitation is 
in any case close to the limit set by our present knowledge of irregularities of 
electron density. 

The forward gain, from Section 2, is 25 for 4»=0-2 and 100 for wy=or!. 
The increased gain at the edges of the beam will not be realized because of the 
variations of 4» with ray direction and because of the necessary frequency 
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bandwidth, nor for the same reasons can interference effects be expected. 
However, a directivity of this order will be of considerable use in the exploration 
of the galactic background, even though the actual value of the gain may be in 
doubt under practical conditions. 

7. Observational methods. (i) Rockets.—It might be possible to observe a 
focusing effect during the short duration of a rocket flight, but the chances of 
success depend on the smoothness of the ionosphere at the place and time of 
launch. Since a rocket traverses the ionosphere rapidly, little will be gained 
by sweeping the frequency of the receiver, and it would be best to use several 
narrow band receivers at separated frequencies, all connected to the same aerial. 
At each frequency there will be two opportunities of observation, during ascent 
and descent. 

(ii) Satelliites.—A satellite, placed in a nearly circular orbit at a height of 
about 400-500 km, would be in a suitable part of the ionosphere for much of 
its orbit. A sweep in frequency over say 2-5 Mc/s, each sweep lasting about 
I minute, would require data to be telemetered at a rate of about 10 numbers 
per second, each with an accuracy of about 5 per cent. This is well within the 
capabilities of typical telemetry systems. Allowing for the difficulties of 
monitoring and recording the telemetry transmission, and for the incidence of 
ionospheric irregularities over as much as half the orbit, there should still be 
about 20 good observations of focusing each day, and a few weeks of observation 
should cover most of the sky with a 20° beam. A narrower beam might prove 
possible for some parts of the sky, but it is unlikely that many observations will 
be possible at a width of 10°. 

(iii) Aerial systems.—Bearing in mind the need for simplicity in satellite 
equipment, it is reasonable to specify that a suitable aerial system must provide 
a signal considerably larger than receiver input noise. An electric dipole can 
be extended to an overall length of 100 ft, either by the extension of flexible 
metal strip or by the release of a thin wire under centrifugal force, and such an 
aerial can produce a signal more than ten times receiver noise even when very 
broadly tuned. The impedance of the dipole will change when it is immersed 
in the ionosphere, but this is of no consequence under the conditions required 
for focusing. 

A magnetic loop, using a ferrite rod of several pounds weight, can also be used 
at frequencies of about 3-5 Mc/s to produce a signal several times receiver noise. 
It is, however, now necessary to tune to the aerial with a narrow-band circuit, 
and the receiver will consequently be more complicated. It is also necessary to 
take into account the variation of inductance of the loop with temperature, 
although there is no effect due to the ionosphere itself. 

In either case the measured field will be inversely proportional to the square 
root of the refractive index. Interpretation of the records will therefore require 
a knowledge of yy. A value of yu») should be obtainable from the record of the 
frequency sweep of the receiver, since a sudden change of output should be 
obtained at the frequency where py) =0. 

Conclusions.—It seems likely that ionospheric focusing sufficient to give a 
forward gain of 25, and possibly more, could be observed at many points in the 
orbit of a satellite, between heights of about 350 to 600km. The chief 
limitations are the irregularities in the ionosphere, and the narrowness of the 
receiver bandwidth that must be used. The frequency range is best chosen as 
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high as is allowed by typical values of the penetration frequency, since the 
bandwidth limitation is less severe at higher frequencies. 

The aerial beam would only apply for one of the two characteristic modes of 
propagation, but it would be sufficiently narrow (say 20° across) to delineate 
some of the major features of the galactic radiation. The principal objects of 
interest here are the spectral indices of the various components, especially the 
remarkable spur of radiation near galactic longitude /7=o0, and also the 
measurement of the distribution of ionized hydrogen close to the galactic plane, 
shown by its absorption of the background radiation. It is unlikely that many 
discrete sources will be detected, since the aerial gain available for focusing is 
inadequate, but it is possible that enhanced radiation from the Sun will be 
observed. 

Observations from rockets are necessarily limited, but would serve to test 
the precision expected from ionospheric focusing. 

Acknowledgments.—1 have been greatly helped in the application of the 
simple theory to a more realistic ionospheric model by discussions with 
Dr K. G. Budden and Mr I. N. Capon, who also computed the ray paths for 
Fig. 6. 
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THE ROTATION OF BARRED GALAXIES UNDER 
GRAVITATIONAL FORCES, II 


S. Aarseth 
(Communicated by F. Hoyle) 


(Received 1961 May 4) 


Summary 


Several models of barred spiral galaxies in rotation under purely gravi- 
tational forces are computed. A homology relation is found thereby 
extending the range of the results available for rotating systems. Approximate 
formulae are given for the period of rotation and the bar mass in terms of the 
inner and outer radius and the central mass. 

The problem is modified by taking the nucleus to be an oblate spheriod 
and it is found that the stability of the bar remains approximately the same for 
eccentricities <o-5. It is predicted that the inner bar will develop trailing 
arms. ‘The mass of the barred system NGC 7479 is estimated to lie in the 
region 2°4—2°9 xX 107° ©, 





1. Introduction.—The problem of the rotation of a barred spiral galaxy under 
purely gravitational forces has been dealt with theoretically in a previous paper (1) 
where one typical model was done by hand calculation. Since then a whole 
sequence of models has been worked out with the help of the digital computer 
EDSAC 2. 

In this problem we take a spherical nucleus with mass M and radius R,. The 
bar is chosen to be cylindrical with diameter 2a, and the mass density per unit 
length is p(r). The two diametrically opposite parts of the bar run from + R, to 
+R,. The equation of motion at an axial point 7, in the bar is given by equation 


(3) of (x); ; 
1 p%= ma, + al | ows — (Ce =e; aya 
- | POH + (a ms ayn} 


+ - Way '— ee an} | (1) 


Let G be the constant of gravitation and T the period of rotation, then 











2n* 
GMT* (2) 
As before we take the density function scaled by M to be 
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The 5 coefficients A,, are determined by introducing the requirement of uniform 
rotation at the points 
n(R,—R 
‘Tp = R, + pach wot er 1) 


5 


where N=0, I, 2, 3, 4 


The parameter « and hence the absolute value of the density function p’(r) 
is determined by introducing the same requirement at r,=R,. The solution 
of p’(r) can be checked, either by using the Stirling-Simpson formula for 
approximate integration or by computing the forces at intermediate points. 

2. Numerical results—The numerical results of the computations for the 
parameter « and the total bar mass M, are given in Table I and some of the 
density functions are drawn in Fig. 1. R, is taken to be 3 kpc and two sets of 
models are included, namely a= R,/3 and a= R,/4. 


TABLE I 


Numerical results of computed models with R,=3 kpc. 


a=4R, a= }R, 

R, a X 10? M,/M aR,*/? a xX 10? M,/M aR,*/2 
6 0°6920 0'233 O°1017 0°6799 o'172 00999 
6°5 0°6166 0282 0°1022 
7 0°5535 0°331 o°1025 
7°5 0°5002 0°382 0°1027 
8 0°4548 0°433 0°1029 0°4425 0°333 O* 1001 
8-5 04158 0°486 0°1030 
9 073818 0°539 O'1031 0°3701* o419* 0°0999 
9°5 0°3522 0°593 0°1031 

10 0°3261 0°648 O°1031 O°3I151 0°508 0°0996 


The model marked * has been computed from a homologous model and eqns. (7) and 


(10). 


If the units for M and T are © and 10° years respectively and the length uni 
is kpc we have from equation (2) 
2_ 4°39 x 10¥ 
a (4) 
Taking M=5 x 10'°©,as a tentative value for the central mass we get T in the 
range I-I to 1°7 x 10° years. 
In columns 4 and 7 of Table I «R,3 is evaluated and it is seen that we have 


aR,22 = ¢ (5) 


where ¢ is very nearly constant for one set of models. Also it is easily seen that 
a linear relation for the bar mass is satisfied to a close approximation, viz. 
M,/M=4dR,+e (6) 

where d and e are constants for the same values of aand R,. The actual deviation 
from the relations (5) and (6) is very small in the examined range of R,. 

In Fig. 1 the density distribution in the bar is drawn for the cases R,=6, 8, 
and 1okpc. ‘The curves for intermediate values of R, are of similar form to those 
shown. 








rm 


on 
or 


he 
he 











No. 6, 1961 Barred galaxies under gravitational forces, I] 


AP (Yh, 


0.06}— 











0.01 
! I ! I | I — 
3 4 5 6 7 8 9 i 
(a) —>" inkpe 
rN 








(b) —>r inkpc 


Fic. 1 


(a) Scaled density distribution in the bar for 3 models with a=}4R,. 
(b) Scaled density distribution in the bar for 3 models wtih a=}R,. 


3. Homology relations.—Two systems are said to be homologous if they have 
the same mass and their size is related by a linear scaling factork. Dimensionally 
we have that GMT? oc L* where L is the length unit. Thus for a homologous 
system of period 7” and scale factor k we have T’?=k*T? and hence from (2) 


a’ =ky (7) 
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where «’ is the parameter corresponding to «. Writing k-! R,’ for R, we then 
have from equations (4) and (5) 


e 12 ° 12 
Tr. 4 ee i Ye R82R, 312, (8) 


For the models a=1/3R, we take c=o0-1027. Using k=R,'/R,=R,'/3 and 
dropping the primes the general formula for the rotation period is 





8-2 12 
T? = 3 (RRs). (9) 


In the cases of a=1/4R, we take c=0-0999 and the constant 8-23 in (9) should 
be replaced by 8-46. Hence T is not sensitive to variations in the bar radius. 

Since M, is the same for homologous models equation (6) gives d'R’,=dR, 
and hence d’=k-'d._ By putting k= R,’/3 and dropping the primes the general 
expression for the bar mass becomes 


34R, 
R, 
For the models a=1/3R, we have d=o-101 and e= —0-375 and in the cases 
a=1/4R, the values are 0-082 and — 0-320 respectively. 

In two homologous models the mass density scales as the cube of the scale 
factor k. In the bar, however, we are concerned with a mass per unit length, 
and this scales as k, viz. p’(r')=k-p(r). It is immaterial in the latter equation 
whether the densities are scaled with respect to mass or not, since the mass M is 
the same for both models. From (3) we then have, using unscaled densities ; 
p'(|r’|) = MZA,,'|r'|" = Mk A, |r|". By equating similar coefficients we get 
A,’ =k-“+¥A,. Again writing k= R’,/3 we can drop the primes and write for 
the general density function 


or)=™ > (2) alt (11) 


The coefficients A,, have been computed for a range of values of R, in the two 
cases a=1/3R, anda=1/4R,. A restricted range of density curves are therefore 
available from equation (11). 

4. Rotational stability.—It now remains to be examined whether the com- 
puted models can retain their shapes for comparatively long time intervals. For 
the bar force at an axial point 7, we have 


Fylty)=_¥ Salty) An (12) 


where the forces f,,(7,,) correspond to terms in p(r). We use the determined values 
of A, and compute f,,(7,) at points intermediate to the fitting points. The bar 
force required for uniform rotation, viz., 


M,/M= 





+e. (10) 


I 


',.4- — 
Pp 2? 
21 p 


has been compared with equation (12) and the deviations of force obtained. 

It was found that the agreement with rigid body rotation is very good every- 
where except near the ends of the bar. In Table II the maximum deviation of 
total force dF/F in the inner and outer bar is given for three representative models. 
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TaBLe II 
Maximum deviation of force in inner and outer bar for R,=3 kpc and a=}R,. 
R, —dF/F Ty dF/F Tp 
6 0°25 per cent 3°24 0°05 per cent 5°8 
7°5 2°2 3°2 orl 72 
9 5°2 3°3 —o'2 8°5 


In the above computations a spherical nucleus has been used. The bar 
mass introduces a force on the matter in the central body which is directed in the 
outward sense parallel to the bar. Hence a more realistic choice would be to 
use a slightly flattened nucleus. ‘Taking an oblate spheroid of constant density 
with major axis R, and eccentricity e the attractive force at an exterior point r on 
the major axis is (2); 

F(r, e)= iz" 1P—sin Bos] (13) 


where 


The model R,=3 kpc, a=1/3R,, and R,=3R,was chosen and density functions 
were computed for a range of e-values. In Fig. 2 p(r) is only drawn for two 
cases; viz. e=o and e=o's since the other curves are very similar. 

It is found that the density function and the parameter « are raised with 
increasing eccentricity, and contrary to expectations the stability in the inner bar 
is not improved. Between e=o and e=o-5 the maximum deviation of the 
required force remains approximately the same and the stability decreases 
gradually above e=o-s. 


A Pls 
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Fic. 2.—Scaled density distribution in the bar for different eccentricity of the nucleus. For e=o, 
a=0°003818, My=0'539M; and, for e=o'5, «=0°004030, M,=0'580M. 
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Equation (5) of (1) gives the time deviation from uniform rotation dT after 

n periods when the relative force deviation is dF/F; 
- af ‘ a . (14) 
If we take dT =o0-1T as the stability criterion we get from Table II nx 400, 200, 
and 100 respectively for the outer bar. The corresponding times are in the 
range 4°5 to 1:5 x 10” years for M=5 x 10%. 

The variation of the lack of force in the inner bar shows that the shortest bar 
is also the most stable. This lack of force will lead to trailing of matter just 
outside the nucleus. Indeed, on some photographs the nucleus seems to be 
wrapped in two symmetrical arms, and their presence has been a puzzling feature 
of galactic structure (3). 

5. Discussion.—The existence of an appreciable number of galaxies with a 
straight bar suggests that the bar is a laterally stable configuration. The trans- 
verse velocity of escape from the bar can be estimated as follows. The potential 
due to a cylinder of length 2/ and with mass per length p is 


adTx 


Vm —2Gpln— (15) 
at a point in the middle when the axial distance is a and a<l. For the velocity 
of escape v we then have; 3v?xV. If we take a=1kpc and 2/=6kpc, this 
gives v~25okmsec-! when M=5x10%©. For comparison, the orbital 
velocity at say 6 kpc is 240 kmsec~. The relative velocity of escape from the 
bar is thus higher than one might have suspected and much higher than the 
transverse components that are likely to exist in the bar. Hence the stars in a 
bar will remain in the bar. 

In none of the computed models do we have that p(|R,|) > 0 which is a require- 
ment of physical continuity. A small mass element situated immediately 
outside the end of the bar will be subject to a smaller attractive force than if it 
were at the end point. Given a circular velocity with period 7 the matter in 
question would lag behind the bar in the subsequent rotation. Thus a small 
amount of matter added to the ends of the bar in order to make p(|R,+dr|)=o0 
would lead to trailing arms without affecting much the conditions in the bar itself. 

Barred galaxies are observed to have trailing arms of matter starting from the 
sudden dispapearance of the straight bar and in some cases they extend to form 
a more of less complete ring. Furthermore, it is also observed that the arms 
emerge from the bar at an angle of nearly 7/2. The composition of the trailing 
arms is not known but analogy with normal spirals suggests a predominance of 
gas and dust whereas the main contents of the bar may well be in the form of stars. 
At the ends of the bar the gas and dust may be squeezed out and trailing arms 
would then develop. The influence of such trailing matter would be to start 
breaking down the bar at its ends. 

6. Application to NGC 7479.—Since this paper was written some obser- 
vations of the barred galaxy NGC 7479 have been published (4). It is reported 
that a linear relation exists between rotational velocity and central distance 
out to the point where the bar turns into trailing arms. The total length of the 
bar is 18-5 kpc for an assumed Hubble’s constant of 75kmsec~! per Mpc. The 
reduced circular velocity at the ends of the bar is found to be 220-°8kmsec~', 
making the period of rotation +257 x 10° years. With T and R, known we can 
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use equations (g) and (10) to determine the total mass of the system provided 
that an estimate can be given for R, and a. 

From the photograph and a drawn isophote R,2>1/3R, and aw1/3R,. 

Equation (9) now gives Mx1-9x10%©. The total mass is then 

M+ M,%2°9 x 106, 
neglecting the mass and influence of the trailing arms. Burbidge et al. got a 
total mass of 2:20 x 10!°© on the simple model of an oblate spheroid of constant 
density rotating about the minor axis. 

The value of the total mass is only changed slightly if we use a different value 
for a or take the nucleus to be an oblate spheroid. Furthermore, choosing 
R,=2'5kpce instead, the total mass becomes 2°4x10!Q©; the dependence of 
the mass upon R, is therefore not very sensitive. By writing T7=27R,v— in (9), 
where v is the circular velocity, it is seen that the derived mass depends linearly 
on the adopted distance of the galaxy. 
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Fic. 7.—Map showing the positions of all sources observed. Open circles indicate 
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n circles indicate sources with angular diameters greater than 10’ arc. 
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